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Abstract 



OO ' Sinai's model of diffusion in one-dimension with random local bias is 

studied by a real space renormalization group which yields exact re- 
sults at long times. The effects of an additional small uniform bias 
force are also studied. We obtain analytically the scaling form of the 
r ^ [ distribution of the position x{t) of a particle, the probability of it not 

£^ ' returning to the origin and the distributions of first passage times, in 

an infinite sample as well as in the presence of a boundary and in a 
finite but large sample. We compute the distribution of meeting time 



of two particles in the same environment. We also obtain a detailed an- 
alytic description of the thermally averaged trajectories by computing 
quantities such as the joint distribution of the number of returns and 
of the number of jumps forward. These quantities obey multifractal 
scaling, characterized by generalized persistence exponents 9(g) which 
we compute. In the presence of a small bias, the number of returns to 
the origin becomes finite, characterized by a universal scaling function 
which we obtain. The full statistics of the distribution of successive 
times of return of thermally averaged trajectories is obtained, as well 
as detailed analytical information about correlations between directions 
and times of successive jumps. The two time distribution of the posi- 
tions of a particle, x(t) and x(t') with t > t' , is also computed exactly. 
It is found to exhibit "aging" with several time regimes characterized 
by different behaviors. In the unbiased case, for t — t! ~ t' a with a > 1, 
it exhibits a scaling, with a singularity at coinciding rescaled posi- 
tions x{t) = x(t'). This singularity is a novel feature, and corresponds 
to particles which remain in a renormalized valley. For closer times 
a < 1, the two time diffusion front exhibits a quasi-equilibrium regime 
with a ln(t — t')/lnt' behavior which we compute. The crossover to 
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a t/t' aging form in the presence of a small bias is also obtained an- 
alytically. Rare events corresponding to intermittent splitting of the 
thermal packet between separated wells which dominate some aver- 
aged observables, are also characterized in detail. Connections with 
the Green's function of a one-dimensional Schrodinger problem and 
quantum spin chains are discussed. 
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I. INTRODUCTION 



Studying non-equilibrium dynamics provides a useful route to elucidate the prop- 
erties of systems with quenched disorder. In addition it is very relevant for exper- 
iments, since most such systems form glassy states with ultraslow dynamics and 
usually do not reach full thermal equilibrium within the accessible time scales. This 
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is the case for a variety of experimental systems such spin glasses [0, random field 
systems f2]||, vortex lines in superconductors [f|||, and domain growth in presence 
of quenched disorder. Despite decades of extensive work, there are still a number 
of unresolved issues in the theoretical description of the dynamics of systems with 
quenched disorder. This uncertainty is due, to a large extent, to the lack of physi- 
cally relevant models for which analytical solutions can be obtained, providing clear 
cut answers to well posed questions. The need for such models is all the more acute 
since it is prohibitively difficult to obtain unambiguous answers from numerical sim- 
ulations when the dynamics is ultraslow, especially since the interpretation is often 
blurred by the absence of precise theoretical predictions. Solvable models, where the 
answers are known, should also provide useful testing ground for numerical methods 
by giving clues on the necessary simulation time scales and averaging procedures in 
disordered systems which are often dominated by rare events. 

Some progress has been made in obtaining analytical solutions for the large time 
behavior of mean field type models 0. Although it is still extremely unclear how 
much these mean field results will carry through to short range, finite dimensional 
systems, one outcome of these works has been to demonstrate the existence of 
several possible large time regimes and to attempt to classify them. This provides 
further motivation to study aging dynamics in a larger class of models, in particular 
to study the possible ways of taking the large time t, t' — > oo limits for correlations 
between configurations of the system at a waiting time t' = t w after a quench at 
t = 0, and a later observation time t. 

Other types of approaches, such as droplet descriptions of the statics and the 
non-equilibrium dynamics of disordered systems ||, make use of domain growth ar- 
guments. These approaches emphasize the leading role of thermally activated pro- 
cesses, which should play an important role in short range models, while mean field 
dynamics may be dominated by other type of collective processes ||. The "coarsen- 
ing" of domain structures evolving towards equilibrium has been studied extensively 
in pure models [10| but little is known rigorously for domain growth with quenched 
disorder. Thus these approaches are still to a large extent phenomenological and 
one would like to find models where solid results about aging in the presence of 
activated dynamics can be obtained analytically. A natural hope for that would be 
to study ID models which could be used as testing grounds for more complex D > 1 
cases which have resisted analytic attack. 

A celebrated ID toy model for glassy activated dynamics is the Sinai model, 
which describes the diffusion of a random walker in a ID random static force field — 
equivalent to a random potential which itself has the statistics of a ID random 
walk |pT| . Although this model (with or without a bias) has been much studied, 
the known analytical results |TT|-|T6| usually concern single time and single particle 
quantities and are technically hard to obtain. It is known that this model without 
a bias exhibits non trivial ultraslow logarithmic behavior, as the walker typically 
moves as x ~ (hit) 2 , as well as several dynamical phases with anomalous diffusion 
as the bias is increased from zero. By contrast, there was until now no exact results 
about two time aging dynamics, despite several mostly qualitative and numerical 
studies [0)0 which found interesting aging behavior in this model. In addition, 



the Sinai model has interesting extensions to many interacting particles, and via 
domain walls, to the Glauber dynamics of ID random field Ising ferromagnets and 
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spin glasses in a magnetic field. 

Recently we have proposed an approach, based on a real space renormalization 
group (RSRG) method, which allows us to obtain many exact results for the non 



equilibrium dynamics of several ID disordered systems |18[]. We have shown that 



it applies to the Sinai model as well as to ID disordered spin models and diffusion- 
reaction processes in Sinai's type of energy landscapes. This RSRG method is closely 
related to that used to study disordered quantum spin chains JTS-El]. The crucial 



feature of the RG is coarse-graining the energy landscape in a way that preserves 
the long time dynamics. In Sinai's model the way to implement the RSRG is very 
direct: one decimates iteratively the smallest energy barrier in the system stopping 
when the time to surmount the smallest remaining barrier is of order the time scale 
of interest. Despite its approximate character, the RSRG yields for many quantities 
asymptotically exact results. As in |21] it works because the iterated distribution of 
barriers grows infinitely wide, consistent with fTTfl . 

The aim of the present paper is to show in detail how the RSRG method applies 
to the Sinai model, allowing one to obtain in a simple way a large number of exact 
results. We obtain a host of quantities such as return and first passage probabilities, 
single time correlations as well as two time correlations of the type that are probed 
in aging experiments. Given the long history of Sinai's model, some of the results 
obtained here have been derived previously, by completely different methods. These 
methods include those from probability theory p5| , pl]j26| -p8| , as well as conventional 
methods of the physics of disordered systems, such as the Dyson Schmidt method, 
replica methods, supersymmetry, transfer matrix etc. |[29|| . Despite that, a large 
number of our results are, to our knowledge, novel. Indeed, as we aim to illustrate 
in this paper, the most interesting feature of the RSRG, besides being simple to 
apply, is that it allows one to obtain all these results (new and old) from a single 
method, while other methods usually allow access to only specific types of exact 
results. As we will explain, the only limitations are the ones usually associated 
with any RG method. First, almost by definition, it only addresses and obtains 
exactly the universal quantities, i.e., the ones which are independent of the short 
scale details of the model. Second, it does rely on the global assumption that the 
starting model is within the basin of attraction of the zero temperature fixed point 
studied here, and is thus not "exact from first principles" . This is not a restriction in 
the case of the Sinai model: because rigorous results already exist from probability 
theory, this last assumption can be considered established. 

In two companion papers [30,31], we will detail the applications of the RSRG, 



given as a short account in [[18], to the Glauber dynamics of disordered spin models 
and to diffusion-reaction processes in presence of quenched disorder. These works 
rely heavily on the Sinai model and the present work. Thus we here give a detailed 
presentation of the results for Sinai model. 

An interesting feature of the RSRG is that it demonstrates in a simple and op- 
erational way how the Sinai model is related to other one- dimensional disordered 
models. More formal derivations of such mappings can also be made in some cases 
via free fermion models. For instance, the quantum XX spin chains with disorder 
and the random transverse field Ising chain (RTFIC) are related via Jordan Wigner 
transformations to free fermion models near half filling with disorder in the hopping 
term. This problem is in turn related, via its expression as a random Dirac prob- 
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lem, to a supersymmetric random Schrodinger operator |32| identical to the Fokker 
Planck diffusion operator associated with the Sinai model 0j3^]. Most of these 
relations have been detailed previously in various contexts (see, e.g. for a review 
H1H IH ) ■ These disordered fermion models have been much reinvestigated recently 
as they provide examples of quantum derealization transitions. It may sometimes 
be useful to recast them in terms of the Sinai model where some quantities have 
a straightforward physical interpretation (e.g., the logarithmic Arrhenius diffusion 
over barriers growing as a random walk gives the logarithmic energy dependence of 
the local density of states). The RSRG demonstrates such mappings for the low 
energy (large time) properties in a very direct way, as we will illustrate. Zero drift 
in Sinai's model corresponds to the self dual critical case in the RTFIC and to 
the antiferromagnetic XX chain f20j , [l9|| , while the zero velocity biased phase [|I3| , P 
corresponds to the Griffiths phase of the RTFIC [^TJ and a dimerized XX chain. 



As we show here, magnetization properties correspond to persistence properties in 
Sinai's model. 

Although the idea of studying random diffusion problem via real space deci- 
mation techniques has been used previously, it has been mostly applied to fractal 
or hierarchical landscapes (see e.g. []57fl ) which are designed for such methods. By 
contrast, here the RSRG emerges from the structure of the zero temperature fixed 
point itself, as the natural way to treat diffusion in a statistically translationally 
invariant disordered system, with no ad-hoc assumptions. Interestingly, a similar 
property arises in the problem of the coarsening of the pure ID $ 4 model at zero 
temperature, which can be treated exactly by successive elimination of the smallest 
domains in the system [38], a method reminiscent of the RSRG studied here. Finally 



note that since [T3[] appeared, several new papers have been devoted to the Sinai 
model H|l. 



A. Outline 



The outline of this paper is as follows. Section |TJ contains a pedagogical intro- 
duction to the real space renormalization (RSRG) approach for the Sinai model. It 
terminates with the explicit expressions for the fixed points of the RSRG ( |11 C| ). In 



Section III we compute the averaged single time diffusion front for the symmetric 



Sinai model in ( [111 A[ ) and with a bias in ( Jill JJ| ) . In Section [TV] we study the returns 



to the origin (persistence properties) of the thermally averaged motion, as well as 
the statistics of the jumps, in the symmetric and biased case. In Section [V], we 
study returns to the origin of a single walker, distributions of first passage time 
and of the maximum position as well as the probability of meeeting of two walkers. 
Section [VT| is devoted to the aging properties of the Sinai model and contains a 
general discussion ([VI A|) , calculations of singular parts of the diffusion front ( |VIB 



the full two-time probability distribution ( |V1 D| ), and the analysis of a simpler case 
(|VI D 2|) . The section terminates with the analysis of rare events and calculation 



of the front in the quasi-equilibrium regime ( |V1F[ ) and fluctuations in the single 
time diffusion front. In Section |VII A| the RSRG is studied in a finite size system; 
equilibration properties, first passage times with boundaries and finite size diffusion 



fronts are computed. Finally in section VIII we obtain the Green's function of the 
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associated Schrodinger operators in [VIII B| . Section |IX| contains the conclusions. 
Further technicalities are relegated to various Appendices. 



II. MODELS AND REAL SPACE RENORMALIZATION PROCEDURE 

A. Diffusion models 

Diffusion in one-dimensional random media has been modeled in three ways, 
which usually lead to equivalent classes of behavior in the large time limit. Proba- 
bilists have often studied models discrete in time and space; for instance, a particle 
on points of a one- dimensional lattice, n, which jumps to the right (n + 1) with 
probability p n and to the left with probability 1 — p n . Physicists on the other hand 
have often considered random hopping models, continuous in time but discrete in 
space, described by the master equation: 

7~" — ~{Jn+l,n ~ Jn.n-l) Jn+l,n = W n +i. n P n — W n n+lP-n+1 (1) 

dt 

W n+ i tn and J n +i,n are respectively the transition rate and the current from n to 
n + and W ntn+ i and J n ,n+i = —Jn+i,n from n + 1 to n. Finally fully continuum 
models, with Fokker Planck equation: 

d t P(x,t) = H FP P = d x D(x)(Td x P - F(x)P) (2) 

have also been studied (with D(x) > 0). 

It is useful to distinguish three classes of disorder (within each description) lead- 
ing to different types of generic ID large time behavior (for uncorrelated disorder). 

(i) detailed balance, random diffusion coefficient 

This corresponds to W H;n+ i = W n+ i in = -Dn, n +i in ([[]) or to F(x) = and 
D(x) a random positive function. It is well known that the large time diffusion 
coefficient is D e ^ = (l/W) -1 for uncorrelated disorder and thus that this model 
exhibits asymptotic "normal diffusion" unless the D nt7l+ i have a broad distribution, 
with a tail near the origin, P{D) ~ D~ a (0 < a < 1). 

(ii) random traps 

This corresponds to W n>n+ i = and W n +i, n = 1/t„. Each site is charac- 

terized by a release time, but the exit is with the same probability 1/2 to the left 
or to the right (the jump probability depends only on the starting point). Again 
this model exhibits asymptotic "normal diffusion" unless the release times have an 
algebraically broad distribution. 

(iii) generic case: Sinai model 

In the generic case one can always parametrize the hopping rates as: 

W n , n+1 = T -l e PEn, n+ l e -kWn-U n+1 ) Wn+1)n = T -l e PE n , n+le -^(U n+1 -U n ) (3) 

where j3 = 1/T and T is the temperature. This can be illustrated as in Fig. 1: there 
is a symmetric barrier E n>n+ i between sites n and n + 1, plus an additional potential 
difference. The barrier E n , n+ i gives the average diffusion coefficient (or attempt 
frequency) on the bond. The "forces" on the bonds are f n +i,n — —(U n+ i — U n ) which 
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represent a local bias. In a finite size system (periodic in the U n ) the expressions 
(HD correctly lead to the Gibbs zero current equilibrium measure, e~ l3Un jZ. 

The main case of interest here and studied by Sinai is that of independent random 
forces. The generic case, for uncorrelated disorder and for distributions of f n and 
E ntn+ i with fast enough decay (e.g. faster than exponential) all belong to the class of 
Sinai's model, which is a discrete time model. A similar potential can be introduced 
for this model: U n — Uq — T 1 ^!^ log(pj/(l — Pi))- One can most easily visualize 
these as Arrhenius motion in a random potential U n which itself performs a random 
walk, either symmetric or biased. This motion has been studied extensively and 
it is known that diffusion is logarithmic x ~ In 2 1 in the symmetric case, sublinear 
x ~ t M for a small bias (/i < 1) and with a finite velocity x ~ Vt for // > 1, where /i 
is related to the asymmetry of the force distributions as defined later. 

W n+i , n 
n n+1 
W n , n+ 1 




FIG. 1. (a) General hopping model; (b) interpretation of (a) with a barrier between 
each site and a potential difference (local bias); (c) model studied here. 

For convenience and to be specific in what follows we will mostly study, as our 
basic model, the random hopping model with the choice -E n ,n+i = 0. We will also 
compare with the discrete time model originally studied by Sinai. However, our 
results are much more general and apply to any model within the locally random 
force class (with short range correlations). 

B. Renormalization method 

1. Definitions and RG equations 

As described above, we consider models of diffusion in ID landscapes in which 
walkers perform Arrhenius diffusion in a potential U n (n is a site index). A "force" 
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variable f n = U n — U n+ i is defined on each bond (n,n + l) (indexed as bond n) and 
as in the Sinai model, the f n are independent random variables with distribution 
Q(f)df. The long-time dynamics in such landscapes are primarily determined by the 
large barriers and deep valleys. Thus we need to be able to focus on these aspects 
of the landscape while eliminating as much as possible the effects of the finer scale 
structure. 

We therefore introduce a renormalization procedure, for a given landscape, which 
will allow us, in this way, to study the asymptotic dynamics. We should emphasize 
that we will apply it mainly to the case of forces independent from bond to bond, 
but it can in principle be applied to any ID landscape. The crucial feature which 
is needed for the RG to yield asymptotically exact results, is that the landscapes 
have extremal values of the potential which grow with length scale. This will make 
the distributions of the renormalized barriers broader and broader. In the case of 
the Sinai model, it is possible to follow exactly the RG flow (because the forces 
remain uncorrelated under the RG) and thus to check a posteriori that at large 
scales the distributions of renormalized barriers are indeed very broad. However, 
the procedure is much more general and would also lead to asymptotically exact 
results for correlated landscapes in which the renormalized barriers become higher 
and higher. The difficulty in such correlated cases is to follow the distributions. Of 
course there are ID landscapes for which the RSRG would not give exact results for 
the diffusion behavior: in particular, bounded potentials which have normal diffusive 
behavior. 

The RSRG procedure on a given landscape is implemented as follows. One can 
first group the bonds with the same sign of the force (see Fig. [I] c), and then 
can start, with no loss of generality, from an "antiferromagnetic" landscape (see 
Fig. |2|) with the f' n alternatively positive and negative) but with a distribution of 
bond lengths U. Our starting model is thus defined by f n = (—l) n F n where the 
F n — \U n — U n+ i\ are the useful variables — called here "barriers" — and the two bond 
variables F, I are chosen independently from bond to bond with an initial distribution 
P(F,l). In the presence of a bias one needs two distinct distributions P + (F,l) 
for "descending bonds" and P~(F,l) for "ascending bonds" (opposing the mean 
force), both normalized to unity. Note that the combining together of consecutive 
descending bonds in this way naturally leads to an exponential tail in the distribution 
P~ and likewise in P + . Such exponential tails in barrier distributions will play an 
important role in the physics and in our analysis. 

We are interested in long times when the behavior will be dominated by large 
barriers and it is on these that we must focus. Our RG procedure is conceptually 
simple: in a given energy landscape it consists of iterative decimation of the bond 
with the smallest barrier Y = F min , say F 2 = U 3 — U 2 = T as illustrated in Fig. 0. 
At time scales much longer than exp(F 2 /T), local equilibrium will be established 
between sites 2 and 3 and the rate for the walker to get from 4 to 1 will be essentially 
the same as it would be if sites 2 and 3 did not exist but 1 and 4 were instead 
connected by a bond with barrier 



F = F\ — i*2 + i*3 



(4) 



and length 



I' = h + h + h- 



(5) 
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We thus carry out exactly this replacement. This preserves the zigzag structure (the 
model remains alternating "up-down" ) and the larger scale extrema of the potential 
since the total length and the extrema of U in the segments are exactly preserved. 
Furthermore if the starting distribution is independent forces from bond to bond, 
this remains so under the RG. One then keeps on iteratively eliminating barriers 
r < F < T + dT thereby gradually decreasing the minimum remaining barrier 
height, T. Note that there is no ambiguity in the case of continuous distributions as 
considered here, as one can always neglect the unlikely events when two neighbors 
or next nearest neighbors are within dT. 




FIG. 2. (a) Energy landscape in Sinai model; (b) decimation method: the bond with 
the smallest barrier F m ; n = F2 is eliminated as three bonds are grouped into one (see 
text). 

The above rules for F and I define the RSRG transformation for arbitrary land- 
scapes. In the case of the Sinai landscape where bonds remain statistically indepen- 
dent one can define 

C — F — r 

and introduce Pp{C = F ~ T, I) and Pf(C — F — T,l) which denote the probabil- 
ities that a ± renormalized bond at scale T has a barrier F = T + £ > F and a 
length I, each normalized by J °° d( J °° dlP^iC, I) = 1. One can then explicitly write 
closed RG equations for these two distributions describing their evolution under the 
decimation represented in Fig. |2|: 

POO 

(dr - d c )P±((, I) = Pff (0, .) *, P r ± (., .) » CtI P r ± (., .) - 2P r ± (C, I) / d/'P^ (0, 1') 

Jo 

POO 

+p r ± (C,0 / ^'(^(0,0 + ^(0,0) ( 6 ) 

where *^ denotes a convolution with respect to ( only and *^ with respect to both 
( and / with the variables to be convoluted denoted by dots. The first term on 
the r.h.s. represents the new renormalized bonds, the second the bonds which are 
decimated as neighbors of the smallest barrier and the last comes from keeping the 
distribution normalized. The total number np of bonds in the system evolves as 
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POO 

dmr = -n r / dl' (P r + (0, I') + P r (0, I')) (7) 

Jo 

We need also to introduce the average lengths 

poo poo 

l T = d( dllp±((,l) 

Jo Jo 

of a ± bond, and the total average length Ir = l T + 7 r of a valley that evolves as 

POO 

d r l r = I r / dl' (P r + (0, /') + P r -(0, 1')) (8) 
Jo 

We have of course that n-p ~ 1/Ir- 

The RG equations (|6|) derived here for Sinai's model are identical to those de- 
rived to study the low energy properties of the random transverse field Ising chain 
(RTFIC) in |^TJ (we choose notations and conventions as in [^TJ) using a pertur- 
bative analysis of the effects of the strongest bonds and fields. The reason for this 
is that the two models are in fact formally related, as mentionned in the Introduc- 
tion. At the level of the RSRG equations, the mapping appears in a very simple 
way: the local random fields and the random exchanges Jk in the RTFIC corre- 
spond to the ascending and descending barriers respectively, through the relations 
F2k/T = — \nhk and F 2 k+i/T = — In J^. We can also identify the renormalization 
scale T in both models. For the diffusion model it corresponds to an Arrhenius time 
scale t = toexp(r/T) to go over a barrier F = T, whereas in the quantum model 
it corresponds to the minimal energy scale of the levels which have been eliminated 
Q = Qoe~ r . The duality between J and h in the RTFIC simply corresponds to 
reversing the average force (i.e., x — > —x) in Sinai's model. As will be discussed 
below, the deviation from criticality parameter 25 in the RTFIC corresponds to the 
parameter /i/T in Sinai's model (see |TB|JT^|) which controls the long time properties 



and the various phases and is defined for the original model with unit length bonds 
by 

(exp(-/i/ n /T)) = 1. 

Zero drift corresponds to criticality in the RTFIC [^TJ, while the biased phase with 
zero velocity p!B|JT^] corresponds to the Griffiths phase of the RTFIC [^Tj] as will be 
discussed below. Note however that the physical quantities of interest in the two 
models can be different. 



2. Effective dynamics and validity of the method 

Throughout the paper, we define the "effective dynamics" as the dynamics which 
consists in putting the particle at time t at the bottom of the renormalized valley 
at scale T = Tln(t/to) which contains the starting point at t — (see Fig. |2|). Thus 
in the effective dynamics the particle does not move unless one of the bonds which 
are the sides of the renormalized valley to which it belongs is decimated, in which 
case it jumps to the bottom of the new renormalized bond as in Fig. |^. Here, to is a 
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non-universal microscopic time scale, which throughout the paper we set to unity by 
appropriate redefinition of time units; we can then use interchangably T and Tint. 

Symmetric case 

This effective dynamics is an approximation of the true dynamics. But within 
the RG approach it can be seen that this approximation becomes better and better 
as T = Tint increases since the distribution of barriers Pr{F) becomes broader and 
broader, as is detailed below. Thus the renormalized landscape consists entirely of 
deep valleys separated by high barriers and with high probability the particle will be 
near the bottom of the valley in which it began. Upon rescaling of space as X = x/T 2 
the effective dynamics of the diffusion front becomes exact as T tends to oo as was 
proven in ref . fill . Indeed, the probability that the walker is close — in a precise sense 
that we discuss later — to its position given by the effective dynamics, approaches one 
at long times. This stronger result has also been rigorously established. JL1],|I^,^£| . 

There is clearly a source of error in the approximation of the true dynamics 
by the effective dynamics when two neighboring bonds have barriers T's that are 
within order T of each other. However, the error introduced by assigning a particle 
to one of two almost-equal-depth neighboring valleys rather than splitting its dis- 
tribution between the two valleys will occur more and more rarely at long scales. 
Furthermore, any such error is wiped out by a later decimation which eliminates the 
two valleys in favor of a deeper valley. These errors thus lead only to subdominant 
contributions to the quantities that we will compute — with the exception of tails of 
certain distributions which are dominated by rare configurations of the lansdscape. 
These subdominant corrections can themselves be estimated via the RG. For in- 
stance the rescaled mean square thermal width of a packet ^ (x(t) 2 ) — (x(t)) 2 (with 
overbars denoting averaging over landscapes) tends to for large T which is its value 
in the effective dynamics, but has 1 /T corrections coming from barrier degeneracies, 
estimated in Section ( |V1 H| ). 

Strong differences between the real dynamics and the effective dynamics can 
appear in some quantities, such as the persistence properties studied in Sections 
flTV| , 0). These quantities are usually in some way nonlocal in time and depend 
on the behavior of the system over time. Even in these cases though, as we show 
in Section ( [TV| , |VJ) , it is possible to compute some of these quantities by a proper 
interpretation and examination of the RG procedure. 

Biased case 

In the biased case with a bias 5 > one finds within the RG that the distribu- 
tion of barriers against the drift is no longer infinitely broad. However if 5 is small 
the barriers remain large enough so that the RSRG remains a good approximation. 
Again, this approximation remains exact, in the same sense as above, in the appro- 
priate scaling limit fixed 5T and x/T 2 (corresponding to the critical region of the 
RTFIC). For a fixed 5 one expects that the thermal packet is spread over several 
deep wells, but when 5 — > the contributions of these few additional wells becomes 
subdominant. 

To conclude this Section, we stress that despite its approximate character, our 
RSRG method allows us to obtain exact results for many quantities both for the 
symmetric and the weakly biased Sinai model. 
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3. RG with one boundary: reflecting or absorbing 



We now consider the problem of diffusion in a semi-infinite one- dimensional 
medium defined as x > 0. In practice there are two main types of boundary condi- 
tions for the diffusing particle: (i) reflecting (the current at the boundary is zero) 
and (ii) absorbing (the probability is zero at the boundary). We show in this section 
how both boundary conditions can be treated by adding to the bulk RSRG specific 
rules near the boundary, which we call boundary RSRG. 

Let us start with the zero bias case and a reflecting boundary. This condition 
can be represented by placing a barrier with infinite potential Uq = +00 at x = 
with U\ finite, as is illustrated in Fig. [| 




FIG. 3. Illustration of the RG in presence of a boundary, (al) Reflecting bound- 
ary conditions: the boundary at site x = can be represented by setting Uq = +00. 
(a2) Renormalized landscape, (bl) Absorbing boundary conditions: the boundary at site 
x = can be represented by setting Uq = —00 (b2) Renormalized landscape, with the 
absorbing zone (see text). 

When grouping bonds with the same sign as in the previous section, the first 
bond will always be descending with an infinite barrier F\ = +00 and a length 
l\. The decimation of the landscape then proceeds as in the bulk case except that 
now the first bond is never decimated and when the second bond gets decimated (at 
r = F 2 ) it simply increases the length of the first bond l[ = h+h + h- One can easily 
see that starting from a landscape where bonds are statistically uncorrelated — with 
a distribution Er(l) for the first bond and Pp(F,l) for all the other bonds — they 
remain so under the boundary RSRG. Upon increase in T, the bulk distribution P 
obeys the same RG equation @ while E satisfies: 
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d r E r (l) = -P T (T)E T (l) + Er{.) *i Pr(T, .) *i P r (.) 



(9) 



The case of an absorbing boundary can be treated in the same way since it 
amounts to setting the potential of the site x = to Uq — — oo. This is illustrated 
in Fig. |^. Thus the first bond will always be ascending with an infinite barrier 
Fi = +oo and a length l\ (and thus cannot be decimated). The rules are thus 
the same as above with the same RG equation @ for the distribution Er(l) of the 
length of the first bond. The interpretation is however different: the first bond 
represent an "absorbing zone" such that any particle starting from a point within 
this zone will be absorbed by the boundary before time T = Tint, while the particles 
starting outside this zone are still "alive" (and outside this zone) at T = Tint (with 
probability asymptotically close to one). 

We note at this stage that this equation coincides with the RG equation for the 
endpoint magnetization in the RTFIC; i.e., with the first exchange being Jo — 0. 
Conversely, a reflecting boundary corresponds to the first transverse field being 
h\ — 0. The equivalence is reversed on the other end of the chain 

In the case of a bias, the probability distribution of the first renormalized bond 
Tp (/) (+ when the bond is along the bias, and — when it is against) satisfies: 

poo poo 

d r E±(l) = P?(0,.)* l E±(.)* l d£P±(C,.)-E±{l) dl'P^^V) (10) 

Jo Jo 

which generalizes equation (H) of the zero bias case. 



C. General analysis of the RSRG equations 



In this Section we recall some results from references 0,^,|T^,^| which will be 
used extensively in this paper, about the large T behavior of the solutions of the 
RSRG equations (|| [10]) and discuss them in the context of the Sinai model. 



1. Symmetric model 

We start with the symmetric Sinai model (zero bias, self-dual) and thus the RG 
of a single distribution P + = P . One first defines the large-scale variance a of the 
potential as 

(Ut-Urf « 2<r|Z i _ i |, 

with the distance from i to j. Since both and Ui — Uj are preserved by the 
RG, a is also preserved and determined by the initial model as 2a = f dff 2 Q(f). 
In the remainder of the paper we will absorb a in / and simply study the case 
o~ = l. The units of length are then 1/a. To obtain the full results one must 
change / — > al (and dl — > adl) in the following formulae. The rescaled probability 
Pr{v> ^) = r 3 Pr(??r, Ar 2 ) in terms of the rescaled variables 

V = C/T, X = l/T 2 

satisfies, when Laplace transformed in A — » p: 
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(Td r - (1 + V )d n + 2pd p - 1)P T (V,P) = Pr(0,p)P r (.,p) * v P r (.,p) 



The fixed point solution (£TJ is found to be 

p( v ,p) = a(p)e- Tlb{p) 

VP 



with a(p) 



sinh( v /p) 



Hp) = Vp coth(y^). 



Thus, taking p = 0, one finds the physically natural result that, due to the occurrence 
of long regions which are predominantly up or predominantly down, the coarse- 
grained probability distribution of barriers in Sinai's model is exponential 



, . 9(F-T) (f-v) 
P r (F) ~ — J -e — r 



with a width which grows as (F) ~ T ~ Tint. The total number of bonds satisfies 
(0) and thus decays asymptotically as n r ~ T~ 2 and the average bond length @ 
grows as ~ T 2 . Since T ~ lnt one recovers Sinai's scaling [ITT1 



x ~ In t. 



(13) 



In the following we will need the explicit form of the distribution P(\) 
d V P( V , A) 



cosh(y/p) 



00 / 1 \ i 00 1 

= £ (n + ij^-D-'e-'^i) £ (-l)-(m+i)«-i<^f (15) 



It was shown in |2l| that the convergence towards the fixed point solution 
P*(ij) = e~ v on the critical manifold (i.e., symmetric perturbations) is like 4 with 
eigenvector Pp 1 ^) = (77 — l)e~ v — corresponding simply to a shift in T — plus other 
parts that decay exponentially in T and depend on tails in the initial distributions. 



2. Biased model 

In the case of the biased model one must follow the descending bond distribution 
P + and the ascending bond distribution P_ which are different. Contrary to the 
previous section, it is more convenient in this case to use the unrescaled distributions 
and variables. In terms of the Laplace transforms Pp(p, I) = Jo + °° dle~ pl P^ 1 , I) the 
equation @ reads: 

(dr - 9 c )P ± (C,p) = P^(0,p)P ± (.,p) * c P%,p) + (P ± (0, 0) - PT(0, 0))P ± (C,p) (16) 



As was shown in ||21|| , for large T the distributions P± take the following form, 
in the scaling regime of small 5 and small p with 5T fixed and pT 2 fixed: 
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P±((,p) = U±(p)e 
«r(p) = 



-C«r(P) 



0*(p) 



p + <5 2 coth pVp + 5 2 ] ^5 



sinh [r \Jp + <5 2 ] 



(17) 
(18) 

(19) 



[Note that here we use U rather than the T of pi).] We will also use the evolution 
equations obtained by substituting (|22"D in the RG equation (|16|): 



d r u±(p) -■ 
drU±(p) 



-U+(p)Uf(p) 
-u*(p)U±(p) 



The distributions of barriers alone are: 

25 



p-(0 



p + (0 



,-2T5 



25 



cxp 



-c 



25 



3 2r<5 



exp 



1 - e~ 2rs 
25 



- 1 



e~" - 1 

and the average lengths of the (±)-bonds are respectively given by: 

1 



(20) 
(21) 



(22) 
(23) 



o . 2 (e 5r sinh(5r)-5r) 
2^ 



e~ 5T sinh(5r)) 



(24) 
(25) 



When 5 



-± 

one has I 



|r 2 . and thus the total average length of a valley is: 



r + I 



/ smh(T 5)\ 



(26) 



where n-p is the total number of bonds. 

The convergence towards the solution (|22| ) has been discussed in [^TJ . The above 
solution (^) thus depends on an "integration constant" 25 = u^ip = 0) —u^(p = 0) 
which is determined by the initial condition, and is proportional to the drift. In |^TJ 
it was identified for small 5 as the ratio of the mean to the variance of the original 
distribution Q(f) of the initial independent bonds of unit length 1; i.e. before 
grouping the bonds together: 



/ 



(27) 



It is useful to introduce a parameter [i defined by the unique non zero solution 
of the equation 



dfQ(f)e 



(2* 



This parameter fi has been introduced previously in [^] and is known to determine 
exactly the various phases of the dynamics of the Sinai model with a bias (x ~ 
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for // > 1, x ~ t for /i > 1). Indeed it is also known in random walk theory |43| to 



control the probability of large excursions against the bias. We now show that we 
can interpret these properties within the RG as associated with the exact decimation 
of the landscape @. 

The distributions P (F) of the barriers F > of the zig-zag landscape (see Fig. 
H) obtained by grouping together the consecutive ascending or descending bonds of 
the original discrete model, are related to the original Q(f) distribution through 



OO 



dFe-^PHF) - 9M (29 ) 

where Q+(s) = dfe~ s f Q(f) and Q~(s) = j^dfe^Qif). The difference of 
potential (F + — F-) of the boundaries of a valley of the initial zig-zag potential has 
Laplace transform 

[ <F^KW[ *F-S>-W-) = ^m^L^ ,30, 

where Q(s) = Q + (s) + Q-(- s ) = dfQ(f) e - sf . The definition of//, Eq. @, 
for the original model is thus equivalent for the initial zig-zag landscape to defining 
fi as the unique non-zero solution of the equation 

(e-^ /T )p 0+ (e^- /T ) P - = l. 

But since the renormalized valleys at scale T are constructed from the valleys of the 
initial zig-zag potential and are statistically uncorrelated, this implies that for the 
probability distributions P r (F) of the renormalized barriers at any scale T, 

/oo poo 
dF + e-» F +' T P+{F + ) J dF_e^-/ T P r "(F_) = 1 (31) 



Using the explicit solutions (|22|) for the distributions of barriers, we obtain in 
terms of = u^(p = 0) the following equation for \i 

uj- \ ( a' \ I (32) 



and we thus obtain that the parameter 25 = «p — wjt parametrizing the RG solutions 
( |22| ) indeed corresponds to the parameter fi/T. Note, however, that the expression 
( |27D is only valid for small 5. 

Thus even away from small //, the RSRG allows one to obtain exact information 
on the structure of the landscape, in particular the behavior of the probability of 
large barriers impeding the drift. For large F we have, from ( p2] ) that (for pos- 
itive 5) P~(C) ~ 26e~ 2S< ' implying that the probability of a large barrier, F, is 
~ exp(— fiF/T). As we will see shortly, this controls the anomalous drift exponent 
/i. 
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3. Boundary fixed point solutions 

The RG equation for the distribution of lengths of the boundary bond E T (l) 
defined in section ( [II B 3|) was given in fllpp . In the Laplace variable with respect to 
length the RG equation fllCf ) reads: 

drE±{p) = E±(p){P?{0,p)P±{p) - P?(0,p = 0)) (33) 

For large T using the properties of the fixed point solution (p2|) for P ± and the 
properties fl20|) of the functions U and u, this can be rewritten as: 

d r lnE±(p) = ^li^lM _ u t (0) = ^(^^±(0) _ ln«±(p)) (34) 

M r (p) 

Finally we find: 

^ _ 4(Q) _ ^ m 

r «r(p) sinh((5r)(VpT52coth[r V / pT^] =F5) 

In the symmetric case 5 = we get 

i? ( ) sinhT^ ^ 

ry^coshr^/p 

whose inverse Laplace transform reads in the rescaled variable 
is 

+oo -, +oo 

£(A) = e~ x ^ n+1 ^ = ^= ("l) m e-^ (37) 



.'VTA 

n=— oo " m=— oo 



III. SINGLE TIME DIFFUSION PROPERTIES IN SINAI MODEL 



In this Section we study the diffusion front using the effective dynamics intro- 
duced in Section ( |II B 2|) . 



A. Single time diffusion front for the symmetric model 

We now consider the (single time) diffusion front, i.e., the probability 
Prob(x, t\x , 0) that a particle starting at x at t = be located at x at time t, for 
the symmetric, zero bias case. At large time t the effective (renormalized) dynamics 
corresponds to moving the particle from its starting point Xo to the lower-potential 
end of the renormalized bond at scale T = Tint that contains xq. This is illustrated 
in Fig 0. In a single environment Prob(x, t\xo, 0) is thus localized near the bottom 
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of the bond — i.e., the bottom of a valley — and the rescaled position x/\v?t has a 
delta function shape at large time. 

One can compute averages over environments, or equivalently over initial condi- 
tions xq (with a spatially uniform measure) in a single environment. The average 
diffusion front Prob(xo + x, t\xo, 0) is obtained as follows. The probability that a 
given bond has length I is Pr(l) = J dr]Pr(r],l) and the probability density that xq 
belongs to a renormalized bond of length I at scale V is lPr{l) / f, lPr(l)- Taking into 
account that the distance \x\ between the starting point xq and the bottom of the 
bond is uniformly distributed on [0, 1], one finds after averaging over I: 

1 f°° 

Prob(x,t|0,0)= / dlP r (l) (38) 

ZJ^Pril) J\ x \ 

Using the fixed point solution (|TT| , |T4] ) with T = Tint, we find that the diffusion 
front takes the scaling form: 

Prob( M 0,0) = ^(^) 09) 
with q(X) = -Y t}L e -^\xmn + i? (4Q) 

n=0 

where we have reinserted a. This coincides with the Kesten-Golosov rigorous result 
l|,^7]] for a Brownian potential, as it should [|TTJ since our method gives exact 



results for properties of the rescaled walks x{i)/ In t. 



B. Single time diffusion front for the biased model 

The case of a global bias (/)q > is described by the RG equations (^) with 
P + 7^ P~ . The fixed point (j22D was analyzed in the previous section. It shows that 
at large scales V the barriers impeding the drift have an exponential distribution 
that does not continue to broaden: 

P r -(F) ~ 25e- 25{F ~ r) e(F - r) (41) 

On the other hand, the bonds along the drift become very long with large barriers: 

P+(F) ~ — e~^0(F - T) (42) 

where Fr ~ ^s e2ST ~ Asymptotically in the RG only barriers impeding the 

drift are decimated, since the barriers to go against the bias are very large. The 



distribution Eq. ([41]) is then simply that of potential drops between the impeding 
barriers. One thus recovers the physical picture |TJJ that Sinai's biased diffusion 
renormalizes onto a directed model with traps (ascending bonds) of release times 
r with distribution p(r) ~ 7- _ ( 1+ w. The average length I r of the descending bond 
distribution Eq. fl24]) yields the anomalous diffusion scaling x ~ t M . 



We now compute the average diffusion front Prob(x, t\0, 0) in the case of a small 
average potential drop per unit length 25 > 0. The argument is as in the sym- 
metric case, except that one must distinguish x > from x < 0, which correspond 
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respectively to the starting point being on a descending, (P + ) or ascending (P~ 
renormalized bond at scale T. One thus uses the formula: 



Probst 1 0,0) = =- 



6{x) / dlP+(l) + 9{-x) / dlP T (l) 

J X J — X 



(43) 



This yields, in the scaling limit where T is large while A = l/T 2 and 7 = TS are 
both fixed but arbitrary, the generalization of Eq(p9|): 



Prob(x,t|0,0) 



a 



T 2 (\nt) 2 



q X 



ax 



T 2 (\nt) 2 



,7 = T8\nt 



(44) 



with 

q{X :1 ) 



7 



sinh7 



6{x)lt;I x - 1 



k cosh k — 7 sinh k 



+9(-X)LT-^ x - 1 



k cosh k + 7 sinh n 



(45) 



with k = a/s + 7 2 . From this expression we can compute the moments. One finds: 
1 



(x(t)) 



r ( sinh [47] — 67 cosh [27] + sinh [27]) 



(46) 



85 2 sinh [7] 2 

(x(t) 2 ) = —-. 1 - . 19 (cosh[67] - IO7 sinh[47] + 3 cosh[47] + I87 2 cosh[27] (47) 
-127 sinh[27] + cosh[27] + 2 7 2 - 5) 



Note that (x(t)) p& f^r 3 for small 7 = ST, a form implied by scaling and analyticity 
in 5. 

One can also perform the Laplace inversion. For 7 < 1 let us introduce the roots 
a t (7) ( n = 0> 1> •••) °f the equation: 



«n(7)cotan(a^(7)) 



with 



nn < 0^(7) < (n + l)n 



(4£ 



For 7 > 1, the root Qfo (7) does not exist, but is replaced by the positive root <5o~(7) of 
the equation 0^(7) coth(«Q (7)) = 7. In terms of these roots, the Laplace inversion 
gives 



q(X,j) = 8(X)J2cMe 



-Xs+(i) 



9(-X)J2cn(l)e- 



\X\Snb) 



n=0 



n=0 



where 



4(7)=7 2 + (<£(7))' 
2 



4 (7) 



7 



2(-l 



,n+l/„,± 



«^(7)) 2 e 



2 P =F7 



,sinh 7 ; v / 7 2 + (^(7)) 2 (7 2 + («±(7)) 2 T7) 
except for the term n = in the domain X > and 7 > 1 for which 

4(7>l)=7 2 -(«o + (7)) 2 



4(7 > 1) 



7 



2(« + (7)) 2 e^ 



sinh 7 ; v / 7 2 -(«o + (7)) 2 (7 + K + (7)) 2 -7 2 ) 
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(49) 

(50) 
(51) 

(52) 
(53) 



Note that 4(7) is an analytic function of 7 despite its definition by two domains. 
In the limit of small 7, we recover of course the symmetric case Eq. (^) using 
a t(l -> 0) = (1 + 2n)|. For large 7 > 1, i.e., Tint > 1/5, we have 0^(7) ~ 
7(1 - 2e~ 2 ~< + ...) and thus 

s+( 7 > 1) ~4 7 V 27 (54) 
4(7 » 1) ~ 4 7 V 27 (55) 

whereas all other coefficients in exponentials are much bigger since 4(7) > 7 2 - m 
the regime 7 > 1, the distribution is thus heavily concentrated to the right of the 
origin and reduces to the simple exponential: 



Prob(z, t\Q, 0) « 6(x) exp[-x/x(t)]/x(t) (56) 
with the mean displacement 



x(t) « t 25T /(4<^ 



We can now compare with known results |11| , |14| : for fixed < /i < 1 the variable 
x = x/t^ 1 is distributed with a half-sided Levy probability density L^x^'^dx" 1 ^ 
where L^z) = LT~f^ z e~ c > iS ' J ' . Our asymptotic result ( [56"D reproduces correctly the 
the small \x limit of this Levy front ]TT|,P1 with the correct prefactor Cu- 



IV. MOTION OF THERMAL AVERAGES: RETURNS TO THE ORIGIN 

AND JUMPS 

We now study "recurrence" properties of the Sinai model. One must carefully 
distinguish between the effective dynamics (i.e., the walker jumping between valley 
bottoms) and the real dynamics. In this section we concentrate on the effective 
dynamics. This amounts, as we will see, to studying the fine structure of the motion 
of the (thermal) packet. Asking similar question for a single particule requires a 
study in the presence of absorbing walls and will be discussed in the next section. 
We will also study the zero crossings of the "running average" 

t 

x(r)dT (57) 

which is an approximation to the thermal average. 

While in a single "run" in a given environment the walker typically crosses its 
starting point many times while trapped in a valley, averaging over many runs in 
the same environment yields a (x(t)) which crosses Xq exactly once each time the 
bond on which x$ lies is decimated, since this causes its valley bottom to cross Xq. 

We will first ask what is the fraction Mfc(t) of starting points, xq, for which the 
thermally averaged position (x(t)\x(0) = Xq) has crossed xo exactly k times up to 
time t. Since the effective dynamics consists in putting the particule at the lowest 
point of the decimated bond, the origin and the particule remain in the same bond 
at all times. The probability of crossing the origin — i.e., the starting point — between 
and t exactly k times is thus the fraction of sites which belong to bonds which have 
changed orientation exactly k times between and t. In particular, the probability 
of no return to the origin is M (t) and is equal to the probability that the bond 
containing the origin has never been decimated. 
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A. Number of returns to the origin: symmetric case 



To compute Mk(t) we use two equivalent methods, which we both describe as 
they will be useful in the remainder of the paper. 

First method: Let N-p(k, rj) be the probability that the bond containing the origin 
has a rescaled barrier 77 = (F — T)/T at T and has switched its orientation k times 
up to scale V. It is normalized as J2t=o Io°° dr]N r (k, 77) = 1 and it satisfies : 

(r«9 r - (1 + v)d n ~ 1) N r (k, 77) = 2P r (0)N r (k, .) *„ P r (.) (58) 
-2P r (0)N r (k, rj) + N r (k - 1, 0)P r (.) *„ Pr(-) 

Introducing the generating function N-p(z,r]) = Y^t=o^ T (k,T))z k , we obtain using 
the fixed point solution Pr(v) = e ~ v 

{Td T - (1 + rj)d v - 1) N r (z, 77) = 2N T (z, .) *„ e ~" - 2N r (z, rj) + zN T (z, O)^"" (59) 

We look for a solution of the form 

N r (z,Tj) = Y-^ z \a(z) + b{z)ri)e- r ' (60) 

and find a quadratic equation ($(z) — l)($(z) — 2) = 1 + z so that 

3- y/5 + 4^ 
= g ( ' 

Second method: this consists of associating with each bond a set of auxiliary 
variables m(k) counting, respectively, the number of sites on the bond which have 
changed orientation exactly k times since t = 0. The RG rules for these variables 
upon decimation of bond (2) read (see Fig 0) 

m'{k) = mi (A;) +m 2 (k - 1) +m 3 (k) (62) 
m'(p) =mi(0) +m 3 (0). (63) 

Introducing the generating function 771(2) = X]fclo m (^) ;2:fc one fi n ds, for a fixed z 
the RG rule 

to' (z) — mi(z) + zm-i (z) + 777,3(2) (64) 

A method to analyze such rules is to write the RG equation for the bond joint 
distribution P(j],m). This is hard to solve, however the RG equation for the first 
moment c(r]) = j m mP(i], m) can be solved. Interestingly, the type of combination 



rule ( pi)) under RSRG has been studied in [glj in the context of quantum spin chain 
models. We have recalled that analysis in the Appendix [A], and generalized it to 
the recursion relation ml = ami + bm 2 + cm 3 , which we will extensively need in 
the present problem. In particular, for z = 0, which corresponds to computing 
the probability of no return to the origin Mo(i), the rule is simply m! = mi + 7773. 
Remarkably, this is exactly the the same as the one for the magnetization in the 
RTFIC |]21||. Thus there is an interesting relation between the magnetization of the 



RTFIC and persistence properties in the Sinai model. 
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Since we are interested in the fraction of initial conditions with k crossings, we 
need the ratio ^j^- where the characteristic length Ir grows as T 2 . The results 



z 



from plj and Appendix |A] is that the ratio ^j^- decays as T *( 2 ) with $ 

(3- y/5TIz)/2 in agreement with the first method (1611). 

Results: returns to the origin and multifractality: From the above result we 
can extract several consequences. First, setting z = we directly find that the 
probability that a thermally averaged trajectory does not return to its starting 
point decays, in terms of l(t) ~ (Tint) 2 , as: 

M (t) ~ ~l{t)-~ e with 6 = 3 ~ (65) 
Second, it is natural to introduce the rescaled number of returns to the origin: 

9 = ^ = W^T) (66) 

and to define the generalized persistence exponent 9(g) characterizing the asymptotic 
decay of the probability distribution of g: 

Prob(^) ~ l(t)-~ e ^ (67) 

We now compute 9(g) from the above generating functional (60). By definition: 



rco pco _ pea 

/ dr]N r (z,r]) oclnT / dg z 9lnr T~ 26 ^ = InT dg e -mr(29(g)-ginz) ^ 
Jo Jo Jo 

Since we know that Nr(z) ~ r~®( z ^ we then obtain, using the saddle point method, 
that $(z) = 29(g*(z)) — g*(z) hiz where g*(z) is the solution of 29 '(g*(z)) = Inz. 
Properties of Legendre transforms thus give that, reciprocally, the exponent 9(g) is 
given by 29(g) = §(z*(g)) + g\iaz*(g) where z*(g) is the solution of $'(z*(g)) = 

2_. We find z*(g) = 2g (g + y/ ' g 2 + 5/4^ and thus: 



«*(ff)' 



%) = fln 



4 2 2 V 4 



The exponent 9(g) is a positive convex function : it decays from 9(g = 0) = 3 ~ 4 V ^ (for 
g = we of course recover the value found previously when studying the probability 
of no return to the origin up to T) to 6>(|) = 0, and then grows again for g > 1/3. 
This implies that: 

k 1 

o = : 7= - — r — > - with probability 1 at large time (70) 
y ln(Tlnt) 3 P 3 & K ' 

All of the moments of g will be dominated by the typical behavior; i.e. (g m ) = 
3~ m for all m. The full dependence on g of the 9(g) function describes the tails of the 
probability distribution of the number of returns of (x(t)), i.e., the large deviations. 

Returns to the origin for the running average: 
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For a given walker, H(£) = | f^x(r)dr will typically behave like (x(t)). We 
conjecture that the probability of k = gln(lnt) sign changes of S(i) up to time 
£ decays with the same exponent 9(g) for g < |. For larger g, the behavior is 
dominated by rare valleys with closely spaced and almost degenerate minima on 
opposite sides of the origin which yield extra sign changes in H(t). We now estimate 
the contribution of these rare events. 

We are interested in situations in which the number of zero crossings N(t) of 
y(t) = JqX(t) dr is much greater than those of (x(t)) whose statistics we know. 
The dominant contributions are from configurations of the random potential for 
which the valley bottom in which the origin lies is split in two halves on opposite 
side of the origin for a very long time. In such configurations the valley has two 
minima at x + > and x_ < with a small free energy difference Te, separated by a 
barrier of some height T . The key point is that if the mean rate of change (dy/ dt) 
{^valley ~ (x + + x_e~ e )/ (1 + e~ e ) happens to be very close to zero, y(t) will change 
sign an anomalously large number of times. To estimate the corresponding number 
of crossings, one can consider, crudely, that y(t) performs a biased random walk, 
with steps of order x±e r °^ T and an average drift er °^ T (^valley (since the typical 
time between jumps is e ro//T ). Using well-known results for biased random walks, we 
can estimate that if the particle is trapped in the double valley for a time r ^> e r0//T , 
the typical number of zero crossings of y(t) in that time interval will be: 

iV(r) ~min( — — , Vre- r °/ T ) (71) 

\Wvalley / 

i.e., with r cutting off a quantity inversely proportional to (x) valley - ^ ne 
tribution of w — ( x ) Ya \\ e y/ x ± is constant near zero, with density of order 1/T , 
because of the distribution of e. Thus we can focus on valleys with the smallest r 
since these will produce the largest number of crossings. 

We must now estimate also the probability that such an atypical valley survives 
for a long time r. For that, we need that neither segment on either side of the origin 
be decimated in the RG for time r, which happens with probability 1/ In 2 r. For the 
contribution of these to the distribution of N(t) we thus have, ignoring constants: 

Prob(iVm = N) ~ / dT [ dw5 ( N - + [ dw5(N - ^) 

J T(lnTf |y T -i/a V w J JO 

~NWW e(Vi - N) (72) 

the dominant contribution coming from the first term in the square bracket. The 
first moment of N is (barely) finite, but higher moments grow with time as: 

t 2 



mr ~ w . (73) 

Thus these type of events completely dominate the distribution of the more-than- 
typical N(t) tail. For g = TV/ In mi > 1/3 the distribution is therefore not multi- 
fractal. 
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On the anomalously small N side of the distribution, the type of events which 
might be troublesome appear to be rare enough not to cause Prob(g) for g < 1/3 to 
differ from that using {x(t}) instead of H(t). The result (p9|) for (x(t)) should thus 
hold also for E(t) for g < 1/3. 



B. Distribution of the sequence of returns to the origin: symmetric case 

We now study a more refined quantity concerning the statistics of returns to 
the origin of the thermal average (x(t)) in the Sinai model. It turns out to be 
possible to obtain the full probability distribution of the complete sequence of the 
times r x = Tlnti, ...T^ = Tlnt^, of successive returns to the origin. This is possible 
because of the remarkable property that every time the thermal packet crosses the 
origin, it "loses its memory" of the past. 

We consider the probability -Dr,r' {v) that a bond has barrier rj at T and has had 
its last change of orientation at scale T'. Its evolution equation reads 

(r«9r - (1 + 7i)d v - 1) D TjT ,( V ) = 2P r (0)D r ,r'(-) *„ M-) ~ 2P r (0)D T , T/ ( V ) (74) 

with the initial condition at T = V given by Dr> t r(v) = Pt'(-) *v -^r'(-)- As stated 
above, this initial condition is independent of previous history because at each dec- 
imation the bond is chosen afresh. Since Pr(v) = e_r '> ^ is natural to look for 
a solution of the form D-p^'iv) = (^r,r' + Br t r'V) e ~ v ■ This is found in terms of 
a = r/r" as: 

Arr = r 4x: ( a ~ X ~ - a ~ X+ ) with A ± = (75) 

Brr = j^y- [(A+ - 1)«' A - + (1 - A_)cT A +] (76) 

The probability for a bond to be decimated at T given that its last decimation 
occurred at T' is p(T,T f ) = — <9 r J °° dr}D r r >(ri) = hp(a = r/r') with: 

p(a) = l-±— (a-*--a- x +). (77) 
a A + — A_ 

Thus we have obtained the probability distribution IT t /(t) for the time t of next 
return to the origin (in the effective dynamics) given that the last return was at a 
time t'\ this exhibits "aging behavior" in a = p4 

„ M) dt i / /int\~ A ~ r\nty x+ \ 

u "^ dt = T^x^r[{^) '■ (78) 

For the sequence of successive returns, the picture we obtain is therefore very 
simple: the sequence of scales {T^} at which the successive changes of orientation 
of a given bond occur is a multiplicative Markovian process constructed with the 
simple rule Tk+i = where {ct^} are independent identically distributed random 
variables of probability distribution p(a). As a consequence, = ak-i(Xk-2 • • • ct^Ii 
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is simply the product of random variables, so that we obtain using the central limit 
theorem that 



lim 

k— >oo 



lnT* 
k 



(In a) 



(79) 



and we thus recover that the number k of changes of orientation grows as lnT = 
lnlnt and that the rescaled variable g = ln ^ lnt ^ is equal to 1/3 with probability 1, 
as in (|70|). 



C. Number of jumps up to time t for the effective dynamics 

In this section we study the behavior of the total number of jumps of the thermal 
averaged position (x(t)) at large times. We introduce the number m(n) of starting 
points on a bond such that the effective — i.e., (x(t)) — walker jumps exactly n times 
between and t for the effective dynamics. We will use m to denote various auxiliary 
variables and trust that such local varying usage will not be confusing. The RG rules 
for these auxiliary variables upon decimation of bond 2 read (see Fig. 0) 

m'(n) = mi(n — 1) + m^in — 1) + m^in) (80) 
m'(0) = m 3 (0) (81) 

Introducing the generating function m(z) = 'Y^ =Q '>Ti{'n)z n one finds the RG rule 
m\z) = zm 1 (z) + zm 2 (z) + m 3 (z). We thus find (see Appendix |A|) that the ratio 
m(z)/lr decays as r~*( 2 ) where &(z) is now the solution of the equation: 



U(-zMz)A) -U(-z,l + Q(z),l) (82) 



in terms of the hypergeometric function U(a, b, z). 

Performing the same saddle point analysis as in the previous section we find that 
the rescaled variable G = n/lnT = n/ln(Tlnt) for the number n of jumps up to 
time t, has a multifractal distribution 

Prob(G) ~ l(t)-^ G) (83) 

where the exponent U)(G) is given by Legendre transform as 2uj{G) = $>(z*(G)) + 
G\nz*(G) where z*(G) is the solution of &(z*(G)) = -G/z*(G). Note that for 
G = (which corresponds to z = 0) one has u(G = 0) = $(z = 0) = 2. As in the 
previous section, the asymptotic value G a that G takes with probability one at large 
times is determined by the minimum of u>(G) where u)'(G a ) = which corresponds 
to z*{G a ) = 1. Thus G a = -z*(G a )&(z*(G a )) = -$'(* = 1). Differentiating (jH) 
and using $(z = 1) = we find: 

G * - £/(-!, 0,l)/2 -(/,(-!, 1,1) - 4/3 (84) 
where Ui(a, b, z) = d a U(a, b, z) and ^(a, b, z) = dbU (a, b, z). 
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A similar method can be used to compute the the joint distribution P(G,g) ~ 
j^-9{G,g) Q £ ^ WQ resca j ec i variables G = nj hi(Thit) and g — kj ln(Tlnt) where 
n and k are respectively the total number of jumps and the number of returns to 
the origin, and hence the associated decay exponent 6(G,g). As an example of this 
application, we give the large time limit (valid with probability 1) of the the total 
rescaled number of jumps G g , conditioned on a fixed rescaled number g of returns 
to the origin: 



G 1 + Z{9) 



S-2$(z(g)) 



fM-i, 1 + $(^)), i) + 1 Z [ 9 \ Z ^ Z{9)) U ^ ^)), I) 



where z(g) = 2g (g + ^Jg 2 + § ) and $(z) = ^vpM. Note that for g = 1/3 one 
recovers Gu 3 = 4/3 as expected. 



D. Correlations of the jumps 

In this section we will obtain some information about the statistical properties 
of the sequence of the directions and times of successive jumps. We will define 
a jump forward as a jump in the same direction as the previous one, and a jump 
backward as a jump in the opposite direction than the previous one. Note that a 
jump backward necessarily involves a return to the origin due to the properties of the 
RG procedure. The directions of successive jumps exhibit strong correlations since 
we have found in the previous sections that the total number of jumps behaves as 
n ~ I mini, whereas the number of backward jumps (returns to the origin) behaves 
as k ~ I In In t. Thus in the effective dynamics the walker is substantially more 
likely to jump in the same direction as the previous jump. This is simply because 
the barrier of a bond which has just been created by decimation and the resulting 
combination of three bonds is higher than that of a typical bond at that scale and 
thus it is less likely to be decimated than the other bond encompassing the valley 
in which the walker rests. 

We first compute the stationary distribution of the number of jumps forward 
made since the last return to the origin, i.e., the probabilities {c p } that a walker (at 
a given time) has made p successive jumps forward since its last jump backward. 
This can be obtained by introducing m(p) as the number of initial points on a 
bond such that the walker has jumped exactly p times since the last passage over 
the origin for the effective dynamics, normalized as Y2^Lo m (p) = I where I is the 
length of the bond. The RG rules upon decimation of bond 2 read (see Fig. |2|) 
m'(p) = m>i(p — 1) + ms(p), m'(0) = I2 + 7713(0). The generating function thus has 
the rules m!{z) = zm\(z) +h + m 3 (z) where I2 = Y^^Lo m ^(p) = m 2( z — 1) is the 
length of bond (2). Similar methods as above then yield the generating function 

c(*) = E^lo c p zP '- 

c(z) = ^-M) + t/(i-^o,i) (g6) 

[) 3(f/(-z,0,l) + f/(-z,l,l)) [m) 

in terms of hypergeo metric functions, normalized to c(z = 1) = 1. From (|86|) one 
gets the c p , e.g. c = (1 + U(l, 0, l))/6 = 0.23394..., c x = 0.17492..., c 2 = 0.13356..., 
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C3 = 0.1029..., C4 = 0.07958... etc. Since c(z) has a pole for z m 1.2884, we obtain 
that 

c p ~ exp(-0.25343p) 

for large p. Also, at any given (large) time the average number of jumps forward 
made since the last jump backward is A^ ward = c'(z = 1) = 3.3975649... 

Next, we study the jump time dependence of backward and forward jumps. In 
Appendix |C 1| we compute the conditional probabilities pp r , (respectively PprO t° 
make a forward jump (respectively backward) at T given that the last jump occurred 
at r'. These are scaling functions of the ratio a — — ^7, i-e., p(, ,b r ,dT = pf ,b (a)da 
with: 

P f (a) = -^{5-(a + 2)e-^) (87) 



2a 3 

p-[a 



'" ■ L(5-(a 2 + 2a + 2)e-^). 



2a 3 

Integrating over a we recover, as expected from ([70] ) and (p4]), the total probabilities 
of the next jump begin a forward or a backward jump as pf = dap? (a) = | and 
p b = |. Note that p(a) = p^(«) + p b {oi) gives the total probability that the next 
jump (/ or b) occurs at Y = aV. 

We now study the statistical properties of the full sequence of the times of suc- 
cessive jumps (ri = Tlnti, . . . I\ = Tint/.). Contrary to the sequence of the times 
of backward jumps studied in Section |IV B| which was simply a multiplicative Marko- 



vian process, there are persistent correlations in the full sequence of jumps which 
makes it much harder to analyze. Indexing each sequence by whether each jump is 
forward (/) or backward (b) we need to introduce the following set of conditional 
probabilities: 

pf- f6 (r fc+1 ,r fc ,...r 1 |r ) (89) 

that, given that a backward jump occurred at To, there are exactly k forward jumps 
occurring at times 1^, . . . Ti before the next backward jump occurs at T^i. These 
conditional probabilities are the elementary building blocks of the full measure for 
the sequence of jumps, since once a backward jump occurs, as was noted in Sec- 
tion |1V El, the process starts afresh. Thus the full measure is simply a product of 



the above terms (p9|) . We have computed the first terms of the set of conditional 
probabilities in the Appendix |C 2j . We obtain, for instance, that, given that the 
previous jump was backward, the probability that the next jump is backward is 
^Ei^—l) = 0.298174 and forward is 0.701826. This result is different from above 
where we did not assume that the previous jump was backward. 



E. Number of returns to the origin: biased case 

We now study the returns to the origin of the thermal averaged position (x(t)) in 
the case where a small bias is applied. Then one expects that the number of returns 
is finite, since eventually the packet will leave the vicinity of the origin. However if 
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the bias is small the number of return is large and universal results can be obtained 
in the limit S small, t large with 7 = 5T\nt fixed. 

The method consists again in introducing auxiliary variables m ± {k) to count 
the number of initial points on a (±) renormalized bond at Y = Tint which have 
changed exactly k times orientation up to time t. The RG rules for these variables 
upon decimation of bond (2) are 

m ± (k) = mf (k - 1) +m^(k - 1) +mf(k) (90) 

for k > 1 and m ± (0) = mf(0) for k — 0. Introducing the generating functions 
m ± (z) = ^2 < k =0 m ± (k)z k one finds the RG rule 

m ± (z) = m 1 (z) + zw^(z) + m 3 (z). (91) 

The calculation of the mean value (m (z)) is performed in Appendix |B| and gives: 

(m±(z)) = 8-*W(A s (j) + e^sinh(7)9 7 ^( 7 )) (92) 

with 7 = 5Y, ip(z) = ^(1 + a/5 + 4i) and A z {^f) = K Z Q ^^(cothj) in term of the 
Legendre function Q u (z). The constant K z depends a priori on z in a nonuniversal 
way. From this we obtain the generating function of the probabilities p 7 (fc) = M^if) 
that the averaged position (x(t)) has returned exactly k times to the origin up to 
time t. It is simply given (since initial conditions are uniformly distributed) as the 
generating function of the total number of initial conditions with k returns divided 
by the total length and thus reads 

f>#>* = <™ + (*)> _+("*"(*)> = 5 ^M Z (Y5) (93) 

with the scaling function: 
2K 

M z (j) = s . nh( ^ )2 [Q^)-i(coth( T )) - coth( 7 )Q; ( , ) _ 1 (coth( 7 ))] (94) 

and normalization implies that K z= \ = 1/2. This function has the follow- 
ing asymptotic behaviors. For small 7 it behaves as a power law M. z {^{) ~ 
K z y/Tr2 1 ~^ z ' ) r^^^] + ^(z))^^' 2 which allows one to recover the results of 
Section [IV A| in the limit of a vanishing bias 8 — > 0. For large 7 — ► 00 it goes to a 
constant Ai z (oo) = K z consistent with a finite total number of returns. 

Setting z = 0, we obtain the probability that (x(t)) has not returned to the 
origin up to time t: 

M (t)=p 7 (fc = 0)~5 2l M (7) (95) 

where the exponent 29 = coincides with the exponent (3 of the magnetization 
of the RTFIC [2T|. Note however, that although the scaling function .Mo (7) in the 



particular case of z = 0, which corresponds to the probability of no return, is closely 
related to the scaling function of the magnetization of the RTFIC, it is not identical 
to it. The probability of no returns of the running average of x, S(t), will have the 
same asymptotic behavior as (x(t)). 
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It is interesting to estimate the distribution of the total number of returns. This 
is achieved by studying the limit r — > oo in (93). We obtain that the probabilities 
Mfc(t) = p-y=oo{k), that the thermally averaged position (x(t)) return exactly k times 
to the origin between t = and t = +00 have generating function 



J>7=oo(A0^ = S 2 ~^M z (oo). 



(96) 



fc=0 



It is thus natural to introduce the rescaled variable g 



(-In §) 



and to look for the 



exponent characterizing the behavior for small S. One finds, by an analysis similar 
to that of Section II V A 



Prob(^) ~ S 2d{9) 



(97) 



with the same exponent 6(g) and multifractal behavior as in (|S9|). The same reason- 
ing as in Section [IV A| leads to the result that g is equal to | with probability one 
for small 5 and thus that the total number k of returns to the origin in the presence 
of a small bias is 



\nS\ 



(98) 



F. Distribution of the sequence of returns to the origin: biased case 

As in the symmetric case, it is possible to obtain the full probability distribution 
of the sequence {r 1; T 2 , . . .} of successive returns to the origin. However in the case 
with drift in direction (+), this sequence is finite with probability 1 since there is a 
finite probability that the particle never returns to the origin if it is on the right of 
its starting point. Therefore the probability Hfc(ri . . . Tk) that the particle returns 
exactly (k) times to the origin from t = to t = 00 and that these returns take 
place at scale (ri, . . . , I 1 *;) can be decomposed as the product 

n fc (ri . . . r 2fc ) = P +(oo, r fc )p"(r fc , r fe _ 1 ) P + (r fe _ 1 , r fc _ 2 ) . . . (99) 

where p ± (r,r / ) are the conditional probabilities that the particle returns to the 
origin at T given that the last return to the origin occurred at scale T' in the direction 
(±), and where p + (oo,r') represents the probability that the particle never returns 
to the origin after the last passage to the origin occurred at T' in the (+) direction 

POO 

p+(oo,r') = 1 - / dTp + (T,T'). (100) 
Jv 

We have computed these probabilities in Appendix [D|. They are most naturally 
written in terms of the reduced variables y = cotli7, y' = coth7 ; with 1 < y < y' 
(where 7 = 5T = T5lnt and 7' = 8V = TSlnt') as: 

P ± (T,T')dT = p ± (y,y l )dy = [Q^^P^y') - Q^y^P^y^dy (101) 
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with = v5j an d Pv an d Qu Legendre functions. One finds, as expected, that 

fi dyp~(y,y') = 1. On the other hand, the probability that the thermally averaged 
position (x(t)) never crosses the origin again after having crossed it at V is: 

p+(oc, = 1- rfI > + ( r > r ') = (102) 

Jt> y + 1 

For 7' — > this probability vanishes as: 

P + (oo,r0 = 2^--| 7 ^ (103) 
0rr(0) 



while it goes to 1 for large 7' as: 



e -4 7 ' 



p + (oo,r') = l-_ = l- 5 ^ (104) 

where /i = 2<5T. The factor l/t' 2fl can be understood with a simple argument. Since 
the particle having crossed the origin at t' belongs to a renormalized bond just 
being created, its barrier is distributed not with Pjt(C) (p2|) , but with Pp, cw+ (C) = 
P+(.) * c P+(.) = (u+,) 2 (e- u t' c ~ ^/f2/i e -C/i/(Tt'M). thig ig depleted near the origin 

which is the key point. For a return to occur after t', the new renormalized bond 
has to be decimated in the future, and the dominant contribution comes from the 
times near t' , i.e., we have to compute the probability that the two independent 
barriers of the neighboring bonds each distributed with P^, are bigger than the new 

(u+,) 2 

(+) renormalized bond at T': this probability is simply ^+ + ' 2u - ^2 an d thus behaves 

as l/(4t' 2fJ- ) at large t'. These events are thus responsible for the dominant behavior 
of p + (oo,r') found above. 



V. RETURN TO THE ORIGIN OF A SINGLE WALKER, FIRST 
PASSAGE TIMES AND MEETING TIME OF TWO WALKERS 



A. Probability of no return to the origin for a single walker 

We now compute, in the presence of a small bias, the probability N + (t) (respec- 
tively N~(t)) that a single walker has remained all the time to the right (respectively 
to the left) of its starting point — the bias being by convention to the right. These 
probabilities are found by placing an absorbing boundary at x = as discussed in 
Section [11 B 3| . We note that the probability distributions E^(l) of the length / of 
the absorbing zone satisfying the RG equation (0) has initial condition 5(1 — 1) 
(counting the first infinitely deep bond in Fig. |3| as length 1 by convention) and it is 
the weight of this ^-function part that determines the desired no-return probability. 
For finite T, E^(l) takes the form: 

r+00 

E±(l) = 5(1 - 1) / d(B±(() + regular part (105) 
Jo 

where R^(C) 1S the probability that the first descending bond (in Fig. |3| for the (+) 
case) has never been decimated up to scale T and has barrier F = F + ( at T. The 



32 



total weight of the delta function rp 1 = J +o ° d(R^(() decays to zero as the regular 
part of E$(l) converges towards the fixed point determined in Section |11 C 3| . We 
obtain the probabilities N ± (t) from A^ ± (t) 



r 



± 

r=Tlni 



The evolution equation for H$(C) reads 

(dr - d c )R$(t) = -Pff (O)iT^(C) + P r T (0)^(.) *c Pr (•) (106) 

Setting -Rp(C) = r r^r (C) we obtain <9rln(rp) = —tip = ^rh^-up") where we denote 
tip = Mp(0) in (|19D . Thus rp 1 ~ which yields: 

N+ V ~ T3^f (107) 

^(*)~^T7 ( 108 ) 

In the biased case, one finds that there is a finite probability N + (t) ~ 25 of never 
returning to the origin if the particle starts in the direction of the drift, whereas if 
it starts against the drift, the probability of not returning to the origin up to time 
t decays as N~(t) ~ 25e~ 25r oc t~^. This corresponds to the probability that the 
origin happens to belong to a "trap" impeding the drift of waiting time larger than t, 



as was discussed in [fL4| . We note that the calculation of these persistent probabilities 



in the Sinai model is similar to the calculation of the endpoint magnetization in the 



RTFIC 21 



In the symmetric case S = we thus find that the probability N(t) that a single 
walker has never crossed its starting point x(0) = x between and t decays at large 
time as 

N(t) ~ \ty e with 9 = - (109) 

with l(t) = (Tint) 2 . Note that the persistence exponent obtained here in the 
presence of disorder is different from the result 9 purc = 1 for the pure diffusion 
problem where the probability of no return to the origin up to time t decays simply 
as ~ It is also significantly larger than the persistence exponent 9 = 3 ~ 4 V ^ = 
0.190983.. for thermally averaged trajectories obtained in Section |VA. 



B. First passage times in an infinite sample 

We now compute the distribution of the first passage time T Xo at x = of a 
walker which start at x = Xq at t = 0. The method consists in placing an absorbing 
boundary at x = and studying the probability S xo (V) that this walker has survived 
up to scale T = Tint as illustrated in (||). We use the method of decimation in the 
presence of a boundary discussed in Section |II B 3| . The probability S% : (T) that the 



walker starting at xq is still alive at T in the presence of a (±) drift is equal to the 
probability that the length l\ of the absorbing zone in near the boundary (see Figure 
|^) is smaller than x a ^ T, which is: 

S±(T)= / dhE^(h) (110) 
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in terms of the function E^ih) studied in Section [II C3| . The probability that the 
first passage time T xo is such that T < T In T Xo < Y + dT is equal to the probability 
o"^ (r) that the walker is absorbed between V and T + dT and is simply obtained: 



<(r) = -drS±(T) = -dr / dlE*{l) 



111) 



(a) 



+ 



o 







FIG. 4. The first passage time T xo at x = of a walker starting at xq is obtained 
from the survival probability in the presence of an absorbing boundary at x = 0. 



1. Symmetric case 

In the symmetric case 5 = we can rewrite ( [[ 1 1| ) in terms of the distribution 
E(X) of the rescaled variable A = given by fl37|) and obtain 

^o(r) = -9r^ dXE(X) = ^E (^) . (112) 

The distribution of first passage time T XQ is naturally obtained in terms of the 
rescaled variable w = T In T Xo / y/xo which is a random variable distributed as 



+00 2 +oo 



U> 3 \W 2 J W 3 ^— ' y/TVW 2 



m=—oo 



This distribution has the following asymptotic behavior. It behaves as s(w) ~ 
(4/u> 3 ) exp[— 7r 2 /(4u> 2 )] for small w as the smaller passage times are strongly sup- 
pressed. However it has a broad tail for large w and decays as s(w) ~ J^ w2 ■ In 



particular its first moment diverges: w = Tin T xo / y/x~o = +oo. One can relate this 
tail to the result of the previous Section concerning the probability N(t) that the 
walker never crosses 0. In general one expects that N(t) ~ C(xq)/T for a walker 
starting at Xq, and that C(xq) is non-universal for fixed Xq. Here we find coming 
from the other limit in the scaling regime Xq/T 2 fixed but small, that the behav- 
ior of C(xo) at large xq should be universal as C(xq) ~ \/xq since we find here 
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2. Biased case 



The Laplace transform with respect to Xq of the survival probability reads 

r^^^mry (114) 

Introducing again the scaling variables 7 = ST, X — f$, the Laplace inversion gives 

00 

^(r) = l-^^(7)e-^ W (115) 



n=0 



where s^ij) have been introduced in (j50|) , and where 

C±( 7 ) = f^L.) MWtfl (11 6) 

except for the term n = in the domain (+) and 7 > 1 where 

C + {1 > 1) = (^L_\ 2 ^)) 2e " 7 ( H7) 

017 ) (7 2 -K + (7)) 2 )(7 + («o + (7)) 2 -7 2 ) 1 J 

For 7 ^> 1, we find using (54) that the survival probability S^(r) in the presence 
of a bias towards the absorbing boundary has a simple exponential dependence on 

S- (r > ~ 1 - e -^ 2 -- 21 = 1 - e -*o/i»/(TV) (U8) 

whereas in the case of a bias in the direction away from the absorbing boundary, we 
find 



CO 



( r » I) - 1 - e- 52 V , 2 ^ 2 ffig 2 2 , e-^ 2 . (119) 
°V 5 J ^ (7 2 + n 2 7r 2 )(7 2 + nH 2 + 7) 

In the limit r — > 00, we obtain the probability that a particle reaches the point 
at a distance Xq from its starting point in the far region against the drift 

S 1 - S »i = ^ e " - ' < 120 > 

which coincides, in the limit /i = 25T — > 0, with the exact result in the regime 
< fx < 2 — corresponding to anomalous diffusion — [^H,H which reads (for T = 1): 

,3/2 r(t) 2 1 -3/2 M 2 Xn/4 
;| r(/u) 1-008(^)^0 e 
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C. Distribution of the maximum position 



The above calculation also yields the distribution of the maximum position 
i max (t) = max 0<i / <t x(t') for a particle starting from x = at t = in the presence 
of a + (i.e., along the positive direction) or — bias: 

/•Xm, 

Prob ± (x niax (0 < x m ) = S*JT = Tint) = / dlE± =Tlnt (l) (121) 

Jo 

and thus the boundary probabilities defined in ( |37l , [3~5"l) Ep =Tlnt (x m ) correspond ex- 
actly to the distribution of the maximum of the Sinai walk. 

In the symmetric case we thus recover via the RG a result derived by Golosov 
pTjj . It is given in explicit form in (j37|). In addition we obtain, in the presence of a 
small bias, the explicit form: 

+oo 

Prob ± (x max (0 = x m ) = E± =TUt (x m ) = ^Q ± (7)4(7)e- XmS " ±(7) (122) 

n=0 

with the conventions of the previous section. 



D. Probability that two particules do not meet up to time t 

In this section we compute the distribution of the meeting time Tj, for two 
particules starting respectively, at x = and x = L. The RSRG method is well 
suited to compute this quantity which may be hard to get by other means. We 
call 1 and 2 the two particles, 1 starting from x = at t = and 2 from x = L. 
We compute the probability ^(r) that the two particles have not yet met at time 
t with T = Tint. The probability 1 — Fl(T) that they have met is equal to the 
probability that the segment [0, L] is included in a single renormalized valley at T. 
The distribution Vp(l) of the length / of the valleys is given as Vp(l) = Pjt *i . 
The probability that both and L belong to the same valley at scale T is simply: 

i r+oo 

l-F L {Y) = r dl(l-L)V r (l). (123) 
tr J l 

This leads to the following expression for the Laplace transform with respect to L 

r dLe- pL F L (r) = i^smM = ( * ) 5hih2 ( r ^;) im 

Jo p 2 l r \psmh {T6)J pcos tf ( T ^/pT¥) + 5 2 

In the symmetric case S = 0, we find Fl(T) — f (A = ]4) with the scaling 
function 

+ 0O 

= E 

— oo 



+oo 



ix'Y] ( 2x + — 



-X7r 2 (n+|) 5 



' l_ e -^ 2 (™+l) 2 \ e - X7r2 ( n+1 z. 
3 ■ —, 2- 



7T 



n 4 



(n+|) 



(125) 
(126) 
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The probability density Hl(T) that the two particles meet between T and T + dT 
is H L (T) = —drF L (T). We thus obtain that the meeting time is T L = exp[w\/L/T] 
where w = Tl/ \fh is a random variable distributed as 

m =' f (l\ 2 £ (* + 1 ,) .-swr (127) 

Note that in the effective dynamics, once the two particle meet they remain 
together at all later times. In the real dynamics, rare events as explored in Section 
|V1H| , can cause them to again split for a limited amount of time in distinct wells 
separated by a distance of order ln 2 t, with a probability of order I /Int. 



VI. TWO TIME DIFFUSION FRONT IN THE SINAI MODEL AND 

AGING PROPERTIES 



In this section we will study the two-time quantity 



P(x,t,x',t') = Probst, x't' 1 00), 

i.e., the probability, over the ensemble of random landscapes and thermal fluctu- 
ations, that a particle starting at x = at t = 0, be successively at x' at t' and 
at x at t. Note that it is normalized as J dx'dxP(x,t,x' ,t') = 1 and thus it is 
different — due to the averaging over the landscape — than the conditional probabil- 
ity Prob(xt\x't' , 00) that the particle be at x at t, knowing that it was at x at t 
and at x = at t = 0. 

The average of the two-time probability contains a lot of information about the 
dynamics after letting the system evolve from t = to t = t' = t w , i.e., the aging 
dynamics. We study P(x, t, x', t') in the limit where both t and t' are large. There 
are several time regimes, according to the precise way that the double limit t', t — > oo 
is taken, and we obtain analytic expressions for the scaling form of P(x,t,x' ,t') in 
each of these regimes. We also study 

Q(y, t, t') = J dx'P{x' + y, t, x', t') , 

i.e., the distribution of the displacements y = x — x' between t' and t. Finally, as 



explained below in [VI A| , we simultaneously obtain results for a "two-temperature' 
evolution. 



Some properties of the quantity P(x, t, x', t') were investigated previously in [pj 



by a numerical simulation and qualitative arguments. Here we obtain instead de- 
tailed exact results for this quantity. Whenever they can be compared these results 
are found in agreement with the conclusions of ||16|| . 

Before presenting the analytical results, let us first give a discussion of the various 
regimes studied in the following sections. 
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A. Discussion of the various regimes 



One can distinguish two main regimes for large t and t', which we discuss in the 
symmetric case. 

(i) "scaling regime": t — t' ~ t ~ t' a , with a > 1: 

This first regime is t ~ t /a with a fixed a = > 1. This regime was called the 
"diffusion regime" in [IB]. In this regime, typically the bond containing the origin 
can been decimated between T' — Tin t' and T = Tint and thus motion can occur. 
P(x, t, x', t'), obtained below by iterating the RG from V to T, takes a scaling form 
in the rescaled position variables. We thus define 

X- x 

and since there are two possible choices for the rescaled x' , we define 

x>- x ' 

and 

Choosing to scale x and x' by the same factor T 2 , the scaling form P a (X,X') for 
the two-time probability distribution reads: 

P(x,t,x',t') ~ — i - A P n _mt [X= — - xo ,X' = — ^ - 1 . (128) 
v ' ' ' ; (Tint) 4 a "hTF \ (Tint) 2 ' (Tint) 2 / v ; 

This diffusion front simplifies in the two limits lnt pa Int' (a = 1) and lnt 3> lnt' 
(a — > +oo). First for r = V one must have: 

P a=1 (X, X') = q(X') 5(X - X') (129) 

where q(X') is the Kesten distribution ( P5] ) obtained previously (note that X' = X' 
for a = 1). An interesting feature is that the delta function component of the 
two-time diffusion front at x — x 1 persists even for lnt > lnt': 

P a (X,X') = D a (X') 5(X - X') + P a {X,X') (130) 

where P is a smooth function of its arguments. This property was suggested in [[L6J . 
Here we find that it arises naturally in the RSRG, since there is a finite probability 
that the bond which contains the origin (starting point) will have its lowest point 
unchanged by the renormalization between V and V (note that the bond can grow 
but only on one side). 

This implies that a finite fraction of particles, D(t, t') = D a = JV/ D a (X') remain at 
the bottom of a valley (their renormalized valley at T') and do not move appreciably 
(i.e., by less than 0(ln t)) between t' and t. 

Finally, for very separated times, i.e., large a, the time evolutions of X at t and 
X' at t' decouple and one recovers again the Kesten distribution (p9|). One has: 
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FIG. 5. Fraction of walkers which do not move between t' and t. 



D, 



,(*') - P™(X, X') - a 2 q(a 2 X') q(X) 



(131) 



In the next section |VI B| we will explicitly compute D a (X'). In section [VI D 
will compute the full smooth part, which is more complicated. 



we 



(ii) " quasi- equilibrium regime": t — t' ~ t' a , a < 1: 
The second regime is for t — t' ~ t' a with fixed a = 



of a is consistent with the previous one in Ref. |L6 



lii(f-f') 

Int' 



< 1. [This definition 



In the second regime, the typical situation is that the thermal packet at t' is well 
equilibrated in a valley with the packet of width of 0(1). In this regime, there is 
typically no motion on scales larger than 0(1) between t' and t as the particle is 
near the bottom of a valley. Motion on larger scales will thus be dominated by rare 
events, which we now analyze. 

First, there is some probability that the valley to which the origin belongs un- 
dergoes a decimation resulting in a jump between t' and t. Although this jump is 
large (the walker will jump to the bottom of a deeper valley a distance of order 
I(r') ~ r' 2 away) the probability that it occurs is of order the probability that 
one of the barriers of the valley at time t' is less than T, which is itself of order 
(r — r')/r' ~ t' l ~ a ~ exp(— (l — a)V) m this regime and thus negligible. The 




FIG. 6. Well with two almost degenerate minima (potential difference of order O(l) 
and a barrier Tq which contribute to the inwell equilibrium dynamics. 

behavior will instead be dominated by rare configurations (but less rare than the 
previous ones) in which the valley at time t' has two almost degenerate minima, 
separated by a barrier r , as represented in Fig. || Jumping between such minima 
persists even for if — > oo. The motion between these minima is equilibrium motion, 
since typically r < T' and the packet are already well equilibrated. In this limit 
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the statistics of the infinitely deep valley potential becomes that of a random walk 
restricted to have — ^valley-min > ^ EUH 

Thus in this second regime t — t' ~ t' a with fixed a = ln [*~, f - < 1 we will 
find that the diffusion front for relative displacements Q(y,t,t') has, in addition 
to a 5-function part for the rescaled variable y = y/\n 2 (t — t') (of weight almost 
equal to 1 which corresponds to the typical valleys), an additional — subdominant 
by j^tj — smooth part: 



fet '^ (rhH))» /a 



(Tln(t-t')) 2 



132) 



where the function f a is universal. This result is obtained in section ([VI F|) . 
(ii) "crossover regime": t ~ t', a = 1: 

Finally, the matching between the regimes (i) and (ii) as t — t' ~ t' is studied in 
Section |\TGj 

Before closing this section, it is interesting to note that by computing the two- 
time diffusion front, we obtain simultaneously the answer to the problem of evolution 
of two independent particles in the same environment, each seeing a thermal noise 
with a different temperature T, for the first particle with trajectory x(t), and T' < T 
for the second with trajectory x'(t). If the two particles start from the same point 
at t = 0, it is clear from the considerations of the effective dynamics that the 
distribution of their respective rescaled positions X = x/(lnt) 2 and X' = x'/(lnt) 2 
will be given by P(x,x',t) ~ P a =T/T'{X, X') with the same scaling function as in 
( J128|) . [Although in the aging problem the thermal noises are identical between t 
and t', this does not make a difference at large times for rescaled quantities.] Note 
that it should be easier to measure the dependence on T rather than Int as the 
latter in practice cannot be varied over a wide range. 



B. Singular part of the two time diffusion front: symmetric case 

1. Probability D a of staying within a well from t' to t 

We start by computing the probability for a particle to stay within the bottom 
of a valley between t! and t: 

/oo 
dX'D a (X'). (133) 
-oo 

We compute the fraction of walkers Dr,r' (C) which are on a bond F = T + £ at 
r and have not moved from T' to T. This means that (i) this bond has not been 
decimated, (ii) one of its neighbors has not been decimated either (the neighbor in 
the same valley, i.e., the right neighbor for a descending bond and left neighbor for 
an ascending one). Thus the bond has been able to grow only on one side. As a 
result D satisfies the RG equation: 

(dr - <9 c )Amv(C) = -Pr(0)£> r ,r'(C) + Pr(0)(P r (.) *c D r ,r>(.) - D r ,r>(0) (134) 
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The second and third terms on the r.h.s come from the allowed decimations of 
the neighbor. Integration over ( shows the total loss of weight from the forbidden 
decimations: (ii) from the first term in the r.h.s. and (i) from the boundary term at 
C = 0. We are interested in large T', and thus we can assume that Pr has reached 
its fixed point value flTT|). Thus one obtains, in the rescaled variables r] = ^ and 



(ad a - (1 + V )d v + i)D a (rj) = / e-KDati-rn). (135) 

Jo 

The initial condition D a= i(i]) corresponds to the probability that a walker be on 
a bond with F = T'(l + rj) at V and, since this probability is proportional to the 
length of the bond, it is obtained from ([TT| ) as 



D a=1 ( v ) = i f°° d\\P(r), A) = i (1 + 2 V ) e~\ (136) 
A Jo 



The solution of (11351) reads 



DM = e-«- 2 (5 + — e- j - e-^^e- (137) 

We note that this barrier distribution is a combination of two exponential factors, 
the expected one exp(— F/T), and the other one exp(— F/V) which represents the 
"memory" from the condition at t'. This will be the typical form for barrier distri- 
butions which we will encounter in all aging calculations. 

We then obtain from (|137|) D(t, t') = D a = J °° dr]D a (r)), i.e., the probability that 



a walker has not moved from t' to t, which takes the remarkably simple universal 
form: 

^-©"(i-^)- (138) 

The behavior of ( p.38| ) for close times a = p4 near 1 is dominated by valleys 
about to be decimated at V. Expanding (|138|) yields D(t,f) ~ 1 — |(a — 1). 
The factor 4/3 is consistent with the most probable number of jumps growing as 
4/31ni found earlier. Let us note that —d t >D(t,t') with D(t,t') given in (|138[) , also 



represents the scaled distribution of the first passage time t' at the bottom of the 
renormalized valley where the particle is at t. This is consistent with the result of 
Golosov p?| . 



2. Weight D a (X') of the delta function component of P a (X, X') 

To compute the full singular part D a (X') in (|130|) we simply have to extend 
the previous calculation keeping track of the length of the bond. Since we are not 
interested in the length at T the length only appears at T' as a parameter in the 
initial condition. The final results reads: 
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D a (X') = a 2 d a (a 2 \X'\) 

One can check that for a = 1, one recovers the Laplace transform of the 

Kesten distribution (|39|) , d\{p') = q(p') = ^7 — f^p 7 }) - ^ so one recovers 

d a (p' = 0) = |A,t' with D(t,t') given in ( |138| ) (with a factor ~ corresponding 
to the total probability restricted to the half axis x' > 0). 

This result ( |139| ) can be explicitly Laplace inverted in the limit of a large a = 
where one can neglect the exponential term, yielding 

o /-+oo 

d a (X') ~ — / dX'X'P(X') = DJ>(X') (140) 

where -D a ~ 5/(3a 2 ) from ( |138| ) is the total weight of particles remaining in their 
wells and ty(X') is the normalized distribution of their positions: 



5 ^ 7r 3 (n + i) 3 \ V 2 

n=— 00 v * \ 

We note that compared to the Kesten distribution q(X' = x' /T' 2 ), $f(X' = x' /T' 2 ) 
has more weight towards the larger values of x'/T' 2 . This is a consequence of the 
fact that the farther the particle goes the more likely it is to be in a deep well where 
it is likely to remain longer without further motion. 

Finally, it is instructive to estimate also the singular part of the averaged condi- 
tional probability Prob(xt\x't' , 00). Using a similar method it is found to be: 

D C a ° nd m = ^ (142) 
where d a (X') is the function defined in (|139|) . 



C. Probability of staying within a well: biased case 

Next we obtain the probability D(t, t') that a walker does not move substantially 
between t and t! (i.e., does not jump to a new valley bottom) in the presence of a 
small bias. This can be computed by extending to the biased case the direct method 
of section |VIB , or from the more general approach presented in the next section. 
Here we only quote the end result: 



D T TI 



sinh 7 



(2 sinh 2 7+I -7'cothY) _ e ~ (7 ~ 7 ' )coth7 ' cosh 7 (cosh 7 



where 



7' 



sinh 7' 



(143) 



7 



5Tlnt and 7' = 5Tlnt'. 



This formula is exact in the small bias scaling limit t,t' —>■ 00, 5 — > 0, with fixed 7 
and 7'. 
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This formula (|143|) is interesting as it exhibits a crossover between two differ- 
ent aging scaling functions corresponding to the symmetric model and the directed 
model, respectively. In the limit 5 — > 0, i.e., 7 — > 0, 7' — > with a fixed ratio 
7/7' = a = Eq. (|138|) is recovered. In the opposite limit in which both 7 ^> 1 
and 7' ^> 1, a nontrivial scaling limit exists when 7 — 7' is kept fixed, i.e., t/t' fixed, 
and the above expression simplifies to: 

D{t, t') ~ e ^ 2(7 ~ 7,) = ^| with L(t) = 7 r ~ t" (144) 

using /i = 25T. This coincides with the small /1 limit of the aging form of the 
directed model (formula (51) of JM) which can be written as D(t,t') = H[{t' /t)^] 

with H[z] = y s *y j f*dy(l — y 1 ^)^^. When /i — > the function H becomes 



exactly H[z] = z, and one recovers (144) 



D. Two-time diffusion front: full analysis 

1. Sketch of the method 



To compute the two time diffusion front Prob(xt, x't'\QQ) in the general biased 
case, we need to introduce quantities associated with bonds which keep track of 
their endpoints and for which a RG equation can be written. We thus define 
fip|, (£, xl, x R ] x'l, x' r ) (resp. Qpp,(£, xl, x r ; x' l , x' r )) as the probability that the 
origin belongs to a descending bond with ends [— x' L ,x' R ] at V and to a descending 
(resp. ascending) bond of barrier ( = F — T and with ends [— xl,x r ] at V. Sim- 
ilar definitions hold with Q + ~ and Q for an ascending bond at V. From these 
quantities one can recover the two-time diffusion front for x' > 0: 

/*oo /*oo /»oo 

Pvob(x, t, x' > 0, f\0, 0) = 9(x) / d( dxL I dx'jQ^'v'iCi x l, x r = x; x' L , x' R = x') (145) 

Jo Jo Jo 

/»oo /»oo /»oo 

+9(-x) d( dx R dx' L Q^,{(,x L = -x,x R ;x' L ,x' R = x') 
Jo Jo Jo 

and similarly for x' < 0: 

^00 ^00 ^00 

Pvdb(x, t, x' < 0, t'|0, 0) = 9(x) / d( dxi / dx' R il^,((, xl, x R = x; x' L = — x' , x' R ) (146) 

Jo Jo Jo 

+9(—x) / d( / dx R / dx' R Q^ l (C,x L = —x,xj i ;x' L = —x',x' R ) 
Jo Jo Jo 

The four RG equations for the four quantities ft^ e r ,((, xl, x r ; x' l , x' r ) and 
flp e r ,((,XL, xr; x' l , x' r ), with e' = ±1 can be written in a compact form: 

(<9 r - <9 C ) n^,(C, x L , x R ; x' L , x' R ) = -2Pff (0)fij^,(£ x L , x R ; x' L , x' R ) 

+ / P^(0,l 2 )P^(.,l 1 )* c fl^,{.,y,x R ,x' L ,x' R )8(x L -(y + l 1 + l 2 )) 

+ / P r T (0, h)P?(., h) * C ^r,r' (•> xl, y; x' L , x' R )5 (x R - (y + l 2 + l 3 )) 

J h>0,h>0,y>0 
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+ 



P r ± (., h) H P r ±(., h)V*? r ,(0, yi, y2; x' L , x' R ) 

ii>0,i 3 >0,yi>0,y 2 >0 

5 (x L - (yi + h)) 5 {x R - (y 2 + h)) 



(147) 



These equations must be solved with the following initial conditions at T = V: 

f2r,' r ,(C, x L , x R ; x' L , x' R ) = 5 e ^5(x L - x' L )5(x R - x' r )lo^(C, x' l , x' r ) (148) 

where Wp,(C, x' L , x' R ) is the probability that the origin belongs at V to a bond (as- 
cending if e' = — 1 and descending if e' = +1) with barrier ( = F — V and of ends 
[—x L , x' R \ at r' : 



u T'{Ci x Li X R 



(x' L +x R )) = ±P^((,x' L + x R ) (149) 



Note that all the primed quantities — those at the earlier time — enter only via 
the initial conditions on the Q'a. These equations ( |147| ), together with the initial 
condition (|148 ) are solved explicitly using Laplace transforms in Appendix |F|. For 
the symmetric case the explicit expression for the Laplace transform of the full 
distribution Prob(xt, x't'\00) with respect to x and x' is given in (|F41| , [F4^ ). In the 
next two subsections we give explicit expressions for some simpler quantities. 



2. Some results for the symmetric case 

We first give the explicit expression for the distribution Q(y,t,t') of relative 
displacements y = x(t) — x(t'), with Q(y,t,t') = J dx'P(x' + y, t, x', t'|00). This 
distribution takes the scaling form: 



Note that we have chosen Y = ^ as the scaling variable here for convenience. 
From Appendix [F]the Laplace transform Q a (p) = J °° dYe~ pY Q a (Y) is found to 

be: 

QM = tanh p { f a) (- 2^ - ^ + (I + 1) coth ^ + sinh ^ ( cosh ^ ~ S ^ 1 ) ) 



H 5— . r r- H 7= cosn VP 

pa A \ smh y/p yjp 



3 VP P-3 



Qpa 2 cosh( v /po;) ycosh^ sinhy/pp— 1 
~ VPCoth ^)(1 + v ^tanh( v ^a)) 
3« 2 (p - 1) 

^ g-Ca-ijv^coth^ / h + 1 ^ sinh^X 

2pa 2 cosh( v /pa) \ cosh ^ sinh ^ y^p ) 

Several properties of this expression can be checked explicitly. First, from normal- 
ization on the half space Y > one has Q a {p — 0) = |. Then, the initial condition 
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at t — t' is Q a =i(p) = \ since Q(y,t,t') = S(y). One can also recover the singular 
part D(t,t')5(y) of the distribution corresponding to walkers which have not moved 
appreciably between t' and t. Indeed one finds Q a (p — > +00) ~ D a /2 where D a is 
given by ( |138[ ). Finally, for very separated times a — > +00 one recovers the Kesten 
distribution (p9|) since the initial condition at V has been forgotten: 

Q«{p=\)^\{l ±-=). (152) 

cr k cosh V k 

From the above result ( |151| ) one obtains the moments of the relative displacements 
which take the general scaling form: 

hit 



(|x(t)-x(tOI»)~(TW)*F h , 
F n [a] =a n (a) + e 1 - a b n (a). (153) 



We give explicitly the form of the second moment: 



, , 61a 4 4a 47 2 409 .„ . 

a2(a) = l8b--T + 60 + ^- 378a^ (154) 

h ^ = -l + i~a + lb (155) 



while the first moment has a simpler expression: 

5 „ 6 221 



< Wt ) - x( t0 |> ~ (T,„0 2 ( ^ - Va + ^ - ^e"<-> ) . (156) 



At large a one recovers Sinai's result (from (p9|)). When t and t' are not too 
separated, i.e., a ~ 1 one finds: 

272 

F 2 [a]~— (a-1) (157) 

and the motion is much slower consistent with the aging of the system. 

Again it is simple to recover the behavior for relatively close times a ~ 1. Ex- 
panding the above (|151|) one finds: 

2Q a (p) ~ 1 - i(a - 1) + (a - l)H(p) + 0((a - l) 2 ) 

cosh y/p sinh y/p \ ^ ^fpj (sinh^) 2 ^ ^ 

The 1 — 4(a — 1)/3, p-independent term represents the probability that the particule 
does not jump and has been discussed in Section VI B . The other term 7i(p) is 
the Laplace transform of the probability that the jump is over a scaled distance 
A = Y, H{\) = J dAi(l + ^)P(0, Ai)P(A - A x ) where the first term corresponds to 
the bond containing the origin being decimated and the second term corresponds to 
the neighboring bond in the same valley being decimated. In Laplace variables this 
gives H(p) = ((—2d p P(r] = 0,p) + P(i] = 0,p))P(p)) which using ( |TTD gives back 
the above result (|158|) . 
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We also give the explicit expression for the normalized dimensionless correlation: 



(x(t)x(t')} _ 72 40 180 2045 



(X*(t))yj(x*(lf)) 



61a 61a 2 427a 3 1281a 4 



/JO M _J6A 

V61a 2 183a 3 61a 4 J v ; 

which decreases from 1 to as a = goes from 1 to +oo. Note that the decay as 
1/at for large a is characteristic of the generic decay of corrections to asymptotics 

as i/r. 



3. Some results for the biased case 



From Appendix |F] one has in principle an exact expression for the Laplace trans- 
formed two-time diffusion front. It is, however, very complicated and thus we give 
here only a few simpler quantities. 

Given that the bias is towards x > and that the starting point is x = 0, the 
probability that the particle is on the side x > both at t' and t has the aging 
behavior: 

P++(r, r') = (6(x(t))6(x(t'))} = 1 / e^e 27 ' - 'H ~ 1) 

16smn 7smn 7 

- 2 7 (e" 27 ' + 2 7 ' - 1) - 2e 27 ' - e^ 7 ' + 4 7 ' 2 + 2 7 ' + 3) (160) 

with 7 = T5 and 7' = T'S, while P (r, V) is obtained by 5 — > —5. The other 
possibilities, P + ~' h (T, T') are obtained from the single time identities P ++ (T, V) + 
P+-(r, T) = P+(T) = Zj/Ir and P ++ (r, T) + P-+(r, r') = P+(r') = Ij,/I r / where 
7 r / are given in ( |24|) 

The explicit expression for the correlation (x(t)x(t')) is given by ( F43|) . 



E. Full two-time aging function in a semi-infinite system 

There is one instructive situation where it is simpler to obtain explicitly the full 
two-time probability distribution Prob(xt, x't'\00), even in the presence of a bias. 
This is the case of Sinai diffusion on a semi-infinite axis (0, +00) with a reflecting 
boundary at x = and possibly some drift in the (+) direction. 

In this case the single time diffusion front is simply Prob^(x't'|00) = E^,(x'), i.e., 
the probability distribution of the length x' of the first renormalized bond near the 
boundary, whose Laplace transform is given in equation (|35|) . Similarly, the two- 
time diffusion front is equal to Prob^(xt, x't'\00) = E^ r ,(x,x'), i.e., the probability 
that the first bond near the boundary has length x' at V and length x at T. Its RG 
equation is given by flTO]) 

poo poo 

drE+ r ,{x,x') = Pf(0, .) * x E+ r ,{.,x') * x / dC'P+iC, .) - E+ T ,{x,x') / dt'Pf(Q, l') (161) 

Jo Jo 
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but with the initial condition 

Ep, r i(x, x) = 8(x — x')E^,{x') 

at r = T'. Note that within the effective dynamics this system has the flavor of a 
directed model, since x(t) — x{t') is always positive. It is thus convenient to define 
the double Laplace transform E^ r ,{p,p') = J +o ° dx'e~ p ' x ' dxe~ p ^ x ~ x '^ E^ r ,(x, x') 
and using the fixed point solution (J22|) for P^r together with the properties ( pCf) the 
above RG equation simplifies into dr In E£ r , (p, p') = dr In U +J > \ . Using the above 

initial condition, we obtain the factorized form E^ r ,(p,p') = E^,(p')Q^ r ,(p) i.e., in 
real space: 



Prob£(xt,x't'|00) = E+(x')Q+ r ,{x - x ) (162) 

where E^,(x') is by ( |122| ) and the distribution of the relative displacements y = 
x(t) — x(t') is 

QtAv) = nrZ, .-"^ Vg+g^tr'^+g) - D (163) 
r>r v y smh(^r) v ^T52coth(rVp+^) -5 

In the symmetric case 5 = 0, Laplace inversion gives the two-time front: 
Q{y, t,t') = ^jrn(n + ~) ^cotan (n (n + ^ ^) e'W^+hr 

n=0 



+ y (%) - || f>mtan (* m fi) J ( 164 ) 



with r = Tint and P = Tint'. 
In the biased case one finds: 



e 7 sinh(7) 



9 +°° 



?1=0 



o +°° 2 2 

+ - — y - m e -^(7' 2 +™v 2 ) 

p/2 75-,^ cotanfvrm^r) — 7' 

m=l v 7 y ' 



(165) 



where the s^(j) are given in fl5"0D and the 6^(7,7') are given in terms of the 0:^(7) 
defined in fl48|) via 

a+(7)cotan ( ^+(7) ) - 7 

a - m) = aM (i66) 

except for the term n = in the domain 7 > 1 where 



(7)coth (^a+(7)j -7 

do(7, 7) = «o + (7) 2 _ + , , 2 V7 2 , ' (167) 

«o (7) -7 2 + 7 
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with given in ([48]). In the limit 7 = T5 = TSlnt > 1 7' = r'5 = Ttflnt' > 1, 
one finds the simple aging form: 

Qt,(«) ~ £w + (1 - £) ^"""p*" (168) 



with /i = 2ST. Equation ( 168|) coincides exactly with the small \x limit of the two- 



time diffusion front in the directed model fj5| . As noted above the model with a 
reflecting wall is in effect asymptotically to a directed model. In the presence of 
a bias, in the long-time limit 7,7' ^> 1 it gives the same results as the full model 
discussed in the previous section, as the influence of the wall vanishes in that limit. 



F. Dynamics within a well 

We now study the dynamics in the time regime t — t' ~ t' a with fixed a < 1. 
As was discussed in [VI A| this is dominated by renormalized valleys at V with two 
degenerate minima U\ and JJ% with JJ\ — U2 of order 0(T) as in Fig. || In a typical 
valley many such degeneracies may exist on small scales y = x — x' of order 1, with 
non-universal statistics, but we are interested in rare valleys where such degeneracies 
exist on scales y ~ T' 2 with barriers T ~ T'; the distribution of these is universal. 

We introduce the probability density i?r(C> h x i To) that a renormalized bond at 
scale T has length /, barrier ( = F — F and has a secondary minimum, degenerate in 
energy with the absolute minimum (i.e., the lower edge of the bond) and separated 
from it by a distance x and a barrier r . The calculation of this quantity is performed 
in Appendix [E[ We find the simple decoupled form: 

1 



with 



Rr(C, I, x, T ) = 9(T - T )0(l - x)P r ((, I - x)^f(x/T 2 ) (169) 

1 



f(X) =^J2 n2 ^ 2 {2Xn 2 n 2 - 3) e~ Xn2 ^ . (170) 



We have written for simplicity R in unrescaled variables, but the expression is of 
course valid only in the scaling regime C~r ~r,a;~/~r 2 (see Appendix for 
details). Note that its total normalization is f Q + °° d( J Q dT f + °° dl f Q dxRr(x, T Q ) ~ 
1/r as expected since it corresponds to a rare event; with T T, the density with 
double minima within T of each other is TR. 

We can now obtain the probability that a Sinai walker will move by y between 
t' and t for t — t'~ t' a with a < 1. We need first the probability Kj>/(y, Tq) that the 
the starting point happens to belong to a renormalized valley at V — T In t' which 
possesses two degenerate minima separated by a distance y and a barrier T < V . 
Taking into account that each of the two bonds forming the valley may be the one 
with the degenerate minima (the probability that both have degenerate minima is 
negligible in the scaling regime of interest) one gets, using ( 170|) 



Kr>(y,T ) = / (h + l2)[Pr'(h)Rr>(k,y,T ) + P r ,(l 2 )R r ,(l u y,T )} 
1 Jh,h 

2(1 + 7^)7^|) (171) 



p 2 ; r 4 Vp2 
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When the starting point belongs to such a valley (characterized by y and r ) the 
walker at T' is well equilibrated (since T < r") and its thermal distribution is, in a 
scaling sense, a sum of two delta function peaks separated by a distance y > with 
weights p = 1/(1 + e~ w l T ) (at the abolute minimum) and 1 — p (at the secondary 
minimum), where w represents the (free) energy difference (of order 0(T)) between 
the two minima. When estimating the distribution of y = \x(t) — x(t')\ one probes 
this equilibrium distribution at t and t' . Denoting T = Tln(t — t') it is clear that y 
can be larger than 0(1) only if T > T in which case it is equal to the separation, y, 
of the two minima with probability 2p(l — p) and small otherwise. Thus to obtain 
the distribution Q(y, t, t') for y in the scaling regime, one must sum over all barriers 
smaller than T = Tln(t—t') (the larger ones contribute only to the already dominant 
delta function part of Q(y,t,t')). Thus, using ( |170| ), we find in the scaling regime 
of fixed a = X ^p- < 1 and y/{Tlnt') 2 : 



Q(y,t,t') = C(T) / dT K r ,(y,T ) 
Jo 



1 + 



V 



(Tint') 



where we have defined: 



T 



(Tln(t-t')) 



(Tln(t-t')) 2 



(172) 
(173) 



n=l v 



1 +°° / 

YV-2 £ ! + 



m=— oo 



2m 2 \ _m£ 



(174) 



Note that the factor C(T) = 2 J +o ° dwe~ w/T /(l + e~ w/T ) 2 = T arises from the fact 
that the distribution of w is constant around w = 0. 

The above result is consistent with previous observations in the case of finite, but 
large, t — t' where moments J dyy k Q(y,t,t') were argued (16| to grow as (Tm(i — 
t')) 2k ~ 1 for k > 1/2 and be bounded for k < 1/2. Here we obtain, in addition, the 
behavior for more separated times with positive (ln(t — t')/ lnt') < 1. 

In the biased case we find similarly (see Appendix [E 1| ) 



6(T - T )6(l - x)P r ± (C, I - x)^f(x/Tl)e 



-xS 2 



r 4 

1 



(175) 



where r(X) is the same function (|170|) as in the symmetric case. From this we 
obtain, as above, the probability K^, (y, T ) that the the starting point happens 
to belong to a renormalized valley at T' = Tint' which possesses two degenerate 
minima separated by a distance y and a barrier r . We find 



r 4 

1 



r 2 

1 



-ys 2 



and thus, integrating over the barriers T < f yields the distribution of displace- 
ments y = \x(t) — x(t')\ in the presence of a small bias. We thereby obtain that 
in the scaling regime where the three scaling variables y/(Tln(t — t')) 2 , 5T In t' and 
a = ln(t — t')/ Int' are held fixed, the distribution of displacements is dominated by 
the rare events with valleys with two degenerate minima and with 
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Q ± (y,t,t')~2 1 



S 2 y 



T 



G 



sinh 2 (5Tlnt') J (Tln(t - t')) 3 V^M* - f)) 



176) 



where, interestingly, G is the same function (|174 ) as in the symmetric case. This is 
because the rare events that dominate are those in which the relevant part of the 
landscape is almost symmetric. 



G. crossover at t ~ t' 



So far we have studied separately the regime t — t' ~ t' a (a > 1) and the 
regime t — t' ~ t' a (a < 1). For completeness let us mention what happens (for the 
symmetric case) in the crossover regime t — t f ~ t' . In this regime, the distribution 
of displacements y = \x(t) — x(t')\ takes the form of a sum of two contributions. 

First, it was found in (|158|) that in the limit a — ► 1 + this distribution is controlled 
by barriers of order F « T. For closer times, r — V ~ 0(T) one must consider more 
precisely the jumping process over the barrier. Asssociating a single relaxation time 
r = e( r ' +<E )/ T with the barrier of height V + e one finds that the contribution of these 
events to the the distribution of y takes the form, (in addition to a piece proportional 
to <%)): 



Q 1 (y,t,t') = f(t,t') 



1 



(Tint 1 ) 



-H{Y 



y 



(Tint 1 ) 



(177) 

The coefficient f(t,t') 
is obtained by noting that the probability that the particle jumps between t' and t 
is {eT 1 ' l T — e~ l / T ) and that the distribution of e is uniform around e = with density 
1/T: 



where the Laplace transform of TC(Y) was defined in 



/(MO 



+oo 



de 



exp(— e e//T ) — exp( 



Tln- 

V 



(178) 



The second contribution comes from the events discussed in the previous subsec- 
tion, corresponding to degenerate wells (Fig. ^j. These are dominant for a < 1 and 
they also give a contribution for a = 1, which must be added to ([177]). The corre- 
sponding contribution Q2^y, t, t') to the distribution of displacements is simply given 
by the limit a — > 1 of equation ( |1 73| ) which corresponds to setting ln(t — t') w lnt'). 

Putting this all together we give the explicit expressions for the second moment 
of y in the various regimes which can be obtained from Q = Qi + Q 2 : 



<(*(*) - x(t')f) 
T(T \n(t-t')f 

T(T\nt') 3 ( 8 1 
(Tlnt') 4 ^ 



45 + 945 



48 fln{t-t') 



8 272, .t.. 

1 hi(-)) 

35 315 y t'" 

hit 



lnt' 



lnt' 

t - t' ~ t' 
t-t' ~ t /Q , a > 1 



t - t' ~ t' a , a < 1 



(179) 
(180) 

(181) 
(182) 



where the function F 2 [a] is given in (|153 ). 
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H. Rare events in the single time-diffusion front 



In this Section we examine further the rare events which produce subdominant 
corrections to the results from the RSRG. As discussed in Section |II B 2| , although 
subdominant, these corrections give the principal contribution to some observables 
when the leading contribution vanishes. This is the case for the thermal width of 
the diffusion front analyzed in |III A|, since this is zero in the effective dynamics. We 



now examine this in the symmetric model. 







(d) 



FIG. 7. Rare events which contribute to to the thermal width of the diffusion front. 
The starting point is indicated by a cross, (a) Valley with two degenerate minima, (b) 
Almost degenerate barriers, (c) Valley just being decimated with a barrier T + e. (d) Rare 
event of higher order with the starting point near the upper edge of a bond. 



The possible rare events which contribute to the splitting of the thermal packet 
are indicated in Fig. |7|. The most important ones, all occurring with probabilities 
of order 1/T, are the following. In case (a) the starting point belongs to a valley 
with two degenerate minima: this is the equilibrium situation already considered in 
the previous Section. In case (b) the splitting is due to the starting point being in 
a valley with two almost degenerate barriers: at scales when the packet overcomes 
the barriers, the packet will split between the two wells located on either sides, an 
intrinsically nonequilibrium phenomenon. Note that if the packet is split at t (as in 
(b)) the probability that it remains split until a later time t decays as (lnt/lnt) 2 . 
In case (c) the walker at T belongs to a valley with a barrier T + e, with e ~ 0(T) 
positive or negative. In this case the thermal packet is already split at T between two 
valleys. There are of course other rare events: for instance the one illustrated as case 
(d), when the starting point is near the upper edge of a bond; this also corresponds 
to an out of equilibrium situation, but it occurs with a smaller probability 0(l/r 2 ). 

Let us estimate in more detail the probabilities Q^\y), Qt\v) an d Qt\v) 
associated respectively with events (a), (b) and (c) in Fig. |7|, that the packet be 
split at T with a fixed separation y between the two parts of the packet. Let us start 
with events (a) and (b) which can be treated similarly. 

Events of type (a): We have already computed in ( |171| ) the probability K-p(y, T ) 
that the origin belongs at T to a renormalized valley having two degenerate minima 
separated by a distance y and a barrier r . Integrating over the barrier r one gets: 

Q { v\y) = £dr K r ( y ,r ) = 2 (1 + ^) ^g (iss) 

where the scaling function G has been introduced in ( |174|) . 

Events of type (b): Here we need to compute the probability Q^\y) that the 
origin belongs to a configuration of type (b) with a distance y between the two 
minima. We first compute the probability R-p(h,yi) that a bond at scale V has a 
length li and two degenerate minima separated by a distance y 1 (|170|) 
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d( J dr R r ((,h, yi ,r ) = P r (l x - yi )-G (^) (184) 



and thus: 



r\ pOO pOO pLi 

Qf{y)=—j d/i jf dl 2 J dy 1 y 1 R r {l 1 ,y 1 )P r (l 2 )6(y-(l 1 +l 2 )) 

y 

1^ dY 1 Y 1 G(Y 1 )(P(.) * P(.)) x= y Yl (185) 
Jo 



2 

f3 



in terms of the scaled P from Eq. (0). 

In either case (a) or (b) the packet is split between two wells and the thermal 
distribution can be written, in a scaling sense, as a sum of two delta function peaks, 
of the form p5(x — X\) + (1 — p)5(x — x 2 ) centered at each minima X\ and x 2 with 
Xi < x 2 , \x 2 — Xi\ = y. In case (a), as before p a = 1/(1 + e~ w / T ) where w is the 
free energy difference between the two minima, while in case (b) a simple estimate 
of the relative escape rates in Fig. |7| also leads to pb = 1/(1 + e~ v ^ T ) where v is now 
the (effective free-) energy difference between the maxima. 

Thus we can estimate the dominant large time behavior of the moments of the 
thermal width coming from the contributions of (a) and (b) which simply add to 
give: 



(\x(t) - (x(t))\*) ia+b) « c k {T) (Tlnt) 2 ^ 1 / dYY k Q^ b \Y) (186) 

where the scaled distribution is 

Q («+6)(y) = 2(1 + Y)G(Y) + 2 / dY 1 Y 1 G{Y 1 )(P(.) * P(.))x=y- Yl (187) 

Jo 

with Y = y/T 2 The coefficients Cfc(T) can be computed using the fact the distribu- 
tions of w and of v have constant density near 0. This gives 

/+OO r> 
dw{e- kw ' T + e- w ' T )/{\ + e- w ' T ) k+1 = -T. 
■oo ^ 

Using the Laplace transformed expression 

-+oo 

dYYQ {a+b) {Y)e- sY 

1 + cosh 2 (v/i) 



2d s v / scoth( v ^) + 



sinh(2 v /i) \sinh(2 v /i) 

We must now consider the events (c) as shown in Fig. [7|. The barrier that the 
particle must overcome to leave the valley is T + e. This corresponds to a single 
relaxation time r = e^ v+e ^ T . Thus the probability that the particle is still in the 
valley at time t = e F//T is simply p c = exp(— t/r) = exp(— e~ € ^ T ). The thermal 
distribution can then again be written, in a scaling sense, as p c 5(x — x±) + (1 — 
p c )5(x — x 2 ) where X\ is the bottom of the valley being decimated and x 2 the bottom 
of the new valley. The distribution of y = \x 2 — xi\ is simply: 
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c\ poo poo poo 

Q ( r c \y)=— / dh / dh / dh(l 1 + l 2 )P r (l 1 )P r (( = 0,l 2 )P T (h)5(y-(l 2 + h)) 
1 Jo Jo Jo 



y 

^ - y 



- / d\ 2 (l + 2X 2 )P( V = 0, A 2 )P(^ - A 2 ) (189) 



o 

The contributions of (c) to the moments thus read: 

(\x{t) - (x(t))|*) (c) « 4(T) (Tint) 2 *- 1 / drr fc Q( c )(r) (190) 

where the scaled distribution reads in Laplace transform: 

p+oa o 1 

y o ^^00= _ £ ^ (mUl ^ + ^-_ ) (Mi, 

Using the fact that the distribution of e is constant near zero, one obtains the 
coefficients dk{T) for k > 1 as: 

/+oo 
rfe [exp(-£;e- £/T )(l - exp(- e - e/T )) + exp(- e - e/T )(l - exp(- e - e / T )) fe ] 
-oo 



k 

ln ( 1 + ^ + E(- 1 ) 1+p ^ ln ( 1+ p) 



P =i 



(192) 



Note that the above argument can be made indentically in the region e < and 
thus we have integrated e from — oo to +oo. 

Our final result for the moments are obtained as the sum of ( |186| ) and ( 190 ) and 
can be computed using the Laplace transforms in ( |188| ) and ( |191| ). Let us give the 
explicit resulting expression for the the lowest moments: 



\x(t) - (x(t))\) w ^(68 + 411n2)T(Tlnt) (193) 



(x(t) 2 ) - (x(t)) 2 w (95 + 681n2)T(Tlnt) 3 (194) 

315 

It would be interesting to measure these quantities in numerical simulations and 
test these predictions. 

Note that the above formula (|186| , |190|) give the leading behavior for moments 
with k > 1/2 which grow with time, while the moments for k < 1/2 are expected to 
be finite and non- universal as in [HJ . This can be compared with the work of Golosov 
[p6[| , who showed the existence of an infinite time limit distribution for y(t) = x(t) — 
(x(t)) and gave an explicit formula in the case of a continuum Brownian potential 
U(x) (which corresponds here to the limit a — > where additional universality 
holds). It is easy to see from ||26|| , as well as from more general arguments, that this 
distribution has a tail l/y 3 ^ 2 at large y. Indeed, from Eqs. ( |174| , |183| , |185| ) we see 
that for t — > oo so that 1 <^ y <^ T 2 , 

gW+W+(c) {y) „ Q (a) (y) ^ (195) 
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Thus moments (y{t) k ) for k < 1/2 should be finite and determined from short scales, 
while moments for k > 1/2 should diverge as t — > oo. Our results include the large 
but finite time behavior and thus go beyond those results of |2^ . 



To conclude, note that the rare events in fl7|) which contribute to the width of the 
thermal packet are also the one which play a dominant role in the aging dynamics 
in the regimes a < 1. We have seen that (a) and (c) are the one which contribute to 
Q(y, t, if) in these regimes. The event (6) does not contribute to Q(y, t, t') (since the 
particle cannot jump back to the degenerate valley) but would have to be considered 
to evaluate P(x, t, x', if) in these regimes as well. 

VII. FINITE SIZE PROPERTIES OF SINAI'S MODEL 

In this Section we apply the RSRG procedure to a finite size system with various 
boundary conditions and obtain exact results for the approach to equilibrium of 
several quantities. 



A. RG for a finite size system 

In order to follow the general measure for a finite size landscape 
one needs to introduce the set of functions representing probabilities 
JV^^; (2, li'i C3> • • • 5 C&-i> h-i', h) in an ensemble of systems of length L with b 
the number of bonds and with barriers Q = Fi — F and lengths U. Note that we 
will not keep track of £1 and (b as these will effectively be ±00 depending on the 
boundary conditions. The normalization condition reads: 

Zl = 2 / ^,^1; C2, h; .-(b-i, h-i; h) = 1 (196) 

Note that under decimation one follows separately b = 2,4,... or b = 1,3,5..., 
depending on the type of boundary conditions studied here either reflecting (R) or 
absorbing (A). Let us write the RG equation for a finite size system, choosing for 
definiteness the case denoted RR below of two reflecting boundaries — as explained 



in section [II B 3| , in this case the first and last bonds have infinite barrier and cannot 



be decimated: 
f 6-1 \ 

k=2 J 

/ 0, 1'; z, 1"; ( 2 , W, ■ ■ ■ ; Cb-i, h-u h) + 

Jz,i+i'+i"=h 
-1 



yj / ^t + l(^ C2, h', ■ ■ ■', Ck-i, h-i] z., I] 0, 1'] (k — I"; (k+i, h+i] ■ ■ ■ ] Cb-i, h-i 

J Z,l + l' + l" = U 



k=2 J Z,l+l'+l"=l h 

rb+2 



+ / K, L (h; (2, Cb-i, h-i; z, I; 0, 1'; 1"). (1 

J z,l+l'+l"=l b 

There exists a quasi- decoupled solution — for Laplace transformed distributions — of 
this equation ( |197| ) (as was also found in the case of the RTFIC J22|), which reads: 
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K >L {h; ( 2 , W,... Cb-i, h-i; k) = E+(h)Pr((2, / 2 )P r + (C3, h)- 

■ ■ ■ P r ((b-2, l b -2)Pf((b-i, l b -i)E r (l b )l r 5 b even ^ 2 



(19£ 



where we have allowed for a bias towards the right, and b is restricted to be even 
since we are dealing with the RR case (see Figure §. In this formula (|9|) P±((, I) 
are the bulk distributions satisfying (f|), /r the average length satisfying (||), and 
the Ep(l) satisfy the semi-infinite boundary RG ([T0|). The integral of the measure 
Zl over all variables in ( |196| ) satisfies: 



w+(0)«-(0) 



1- P-(p)P+(p) u+(p)u-(p) -U+(p)U-(p) p 



(199) 



where we have used ( |22| , |35|), and thus Z L = 1 so the finite size measure is correctly 
normalized. 




FIG. 8. Schematic of a finite size system of fixed length L with reflecting boundaries 
(case RR). 



In the case of two absorbing boundaries AA, the solution of the corresponding 
RG equations is obtained by simply exchanging + and — in ( |198| ) (b remains even). 
In the case RA the solution reads: 

Kjh; C2, k; ■ ■ ■ Cb-i, k-u k) = E+ (h)Pf (C2, / 2 )Pr + (C3, h) (200) 
• • .P^(Cb-2,lb-2)Pf (G-i,l b -i)E+ (l b )l r 5 |l b odd ^ 3 ( 201 ) 

together with the term 6=1, which corresponds to the final state with a single 
(+) bond over the whole system (all particles having been absorbed by the right 
boundary), and has for probability: 
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N^ L = [ dV [ iV* >L (Z l5 C 2 = 0, l 2 ; l 3 ) (202) 
Jo Jh,l2,l 3 

Finally, the solution in the AR case can be obtained by exchanging + and — in 

(HD- 



B. Evolution towards equilibrium in a system with reflecting boundaries 

Let us start by studying the equilibrium and the approach to equilibrium in the 
system of size L with two reflecting boundaries. The equilibrium state corresponds 
to the large T limit of the measure ([1961) . In that limit only the term 6 = 2 ( p8|) 
survives and it corresponds to equilibrium in a single renormalized valley. It will 
be reached, as illustrated in Fig. || when the last bulk bond is decimated at some 
T = T eq = T\nt eq , with a certain sample to sample distribution for the equilibrium 
time t eq , which we now compute. 



1. Distribution of equilibration time 



The probability for a sample to reach equilibrium between Y eq and T eq + dT eq , 
i.e., the probability that the slowest relaxation time t eq be such that T eq < Tlnt eq < 
r eq ~t~ e g, is 



PL (T eq ) = d Teq / Nl q>L (h, l 2 ) = d Teq Ur^O *l (203) 
and using (|35|) , the Laplace transform with respect to the system size L is 



J^dLe- pL p L (T eq ) = d req ^ 



(j9 + 5 2 )coth 2 T eq ^p~T¥ 



-5 2 



For zero bias, one introduces the scaling variable 

Tint 



w 



>'(:(/ 



and finds that it is distributed as: 



+oo 



P(w) 



£ 2 * 2 

n=— oo \ 

+oo 



n 



VI 



i\ i \ 

1 e 



2 J w 2 



V (-l) m (l -2mV)e- mV 



In the presence of a bias one can compute, e.g., the average: 



T In t eq 



2W- rdue-^lnt 1 - 

7T ./ V If 



(204) 



(205) 
(206) 



(207) 



56 



2. Distribution of equilibrium position x 



cq 



The probability that the bottom of the single remaining equilibrated valley is at 
x = can be obtained as 

QlM = N^ L (x eq , l 2 ) = (l r E+(x eq )Er(L - x^))^ , (208) 

which leads simply to: 

Qdxeq) = e + (x eq )e~(L - x eq ) (209) 

where 

e ± (x) = LT^ X {^—) = LT£ s -=±=— (210) 
«r=oo(p) \/p + 5 2 =f 5 

Thus we get e + (x) = e~(x) + 25 with 

e"(x) = -Le-* 52 - f£ / dv 6 — — . (211) 



At small x one has e ± (x) ~ ^= while for large x, e + (x) ~ 25 and e (x) ~ e z<52 . 
Thus in the biased case with L S> 1/5 2 in equilibrium the particle is confined within 
a distance y = L — x eq ~ 1/5 2 near the left boundary distributed as 2<5e~ (?/) . In the 
symmetric case 5 = the equilibrium position is distributed over the whole system 
as: 

Qdx eq ) = 1 (212) 

which has a simple probabilistic interpretation in terms of the landscape random 
walk confined to U(x) > U min = U(x eq ) on both sides of x eq . 

Finally, we obtain the joint distribution of equilibrium position x eq and equilib- 
rium time T eq = Tlnt eq : 

P L (T eq ,x eq ) = 8 Teq (l Teq E + (x eq )E-(L - x eq )) (213) 

where E ± was computed in (|35|j37|). 



C. First passage times 

1. With a reflecting boundary 

Let us compute the probability S XOj l(F) that a walker starting at xo is still 
alive at T in the presence of an absorbing boundary at x = and a reflecting 
boundary at x = L. It can be expressed as an average over the measure AR (|196 , |200|) 



S X o,l(^) — (@( x o — h))- Thus its Laplace transform with respect to L reads 

/ dLe- pL S XOiL (T) = -p?L / dhe- ph E r (h). (214) 
Jo pEriP) Jo 
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In the particular case where the starting point coincides with the reflecting 
boundary (x = L) it is simpler to obtain the first passage time T LL at x = 0. 
In that case, the probability to be absorbed coincides with the probability pz,(r) 
that the last decimation occurs in the AR system at T; from (|202|) : 

p L (T) = d r N^ L = l r E+(.) * L P r "(0, .) * L E+(.). (215) 

In the symmetric case we thus obtain that the scaled first passage time variable 
w = T ln Jj^ is distributed as: 

s(w) = ^ £ (-1)" (n + I) e^Ki) 2 = ^ E (-!) m ( m + I) e~< m+ ^ ■ (216) 



w 3 ^ v ' V 2 7 v/tt ^ V 2 

n=— oo x ' v m=— oo x ' 

We note that the distribution of T^x, was obtained previously by a completely 
different method in [ £Efl . Here we also recover ln^T/^) 9 ) ~ g 2 L, i.e., that the first 
passage time is a strongly fluctuating quantity 



2. With an absorbing boundary 

Let us now consider an absorbing boundary at x = L. We first compute the 
probability po(xo, L) that the walker starting at xq reaches x = before x = L. Since 
the final state of the AA system consists of two absorbing zones associated with each 
boundaries, the first one [0, x eq ] and the second [x eq , L] where x eq is distributed as 
in ( |212| ). Thus, in presence of a bias applied in the direction ± the result reads: 

p±(x ,L)=f dxe T {x)e ± {L-x) (217) 

Jx 

where e=F(x) was computed in ( |211| ). In the symmetric case this gives po(xo,L) = 
^Arccos (2f - 1). 

One can also compute the survival probability S XQ} l(T) of a walker starting at 
Xo in the presence of two absorbing boundaries at x = and x = L. It is obtained 
as an average over the measure (|E5|) of the finite size system AA as 



and thus reads: 



S X0 , L (T) = (^o - h)0{L -k- x )) (218) 



c (r^-TT- 1 tanh(y^+g) tanh(^/p) / 1 

XoM J L1 r+L/r**+*ofl* q^>^p-+-q ^sinh 2 ^) smh^P + ^J ' ^ 

Note that the distribution of times ti ast at which all particles have left a given AA 
sample is identical to the one computed in (|206|) as lnti ast = lnt eg . 
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D. Averaged diffusion front 



We first discuss a system bounded by two reflecting walls at x = and x = L. 
The full averaged diffusion front Proboi(a;, t\xo, 0) for a walker starting at xo at t = 
is computed in Appendix for both biased and symmetric cases. 

In the symmetric case it takes a scaling form Proboi(x, t\xo, 0) = f^qx(X\Xo) as 
a function of the rescaled variables X = x/T 2 , X = Xq/T 2 and the rescaled length 
of the system A = L/T 2 where T = Tint. The Laplace transform q{p,Po,q) = 
jo + °° dX Jo dX Jo dX e-( pX+PoXo+qX kx{X\X ) of the rescaled front is: 

q(p,Po,q)= P P+Po+ g P P + q -Ppo+ q P l 1 _ Epo+g (220) 

Poq(p + Po + q)E p+po+q E q Poq E q Po{P + Po + q) E p+po+q 

where P p — 1/ cosh(y / p) and E p = tanh( v /p)/ v /p. In the limit L — > +oo we can 
obtain the averaged front in a semi-infinite space with a reflecting wall at x = as: 

Km (9 qWo,9)) = P ; +P ° Pp r/ P ° + -4- (221) 
^° Po{p + Po)E p+po Po 

The corresponding formula for AA and RA are given in Appendix 0. 

There is a case where the diffusion front in a finite sample takes a particularly 
simple form. This is when the starting point coincides with the reflecting boundary 
Xo = 0. The calculation of the Appendix simplifies as one then has that x — l\ 
where l\ is the length of the first bond. In the RR (or RA) case: 



Prob L(x, t\x = 0, 0) = £+(x)0r' A (L - x) (222) 



with 



1 

'pEr(p) 



(*) = ^T^TTTT (223) 



{x) = LT^ x ^M- (22 H 

pE r (p) 



In the symmetric case one finds simply: 
-if/ \ rrr-i r coth(r^p) 



[X 



LT;i ;—y" =2^e"^^ (225) 



n=0 
oo 



y r \x) 



^vpsto^S^* 8 "^*- (226) 



VIII. RESULTS FOR FOKKER PLANCK AND ASSOCIATED 
SCHRODINGER OPERATOR 

It is also interesting to obtain results, via the RSRG, for the random Schrodinger 
operator associated with the Sinai diffusion problem. We first recall the connection 
between these two problems. In this Section we set T = 1. 
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A. Prom Fokker Planck to Schrodinger operator 

In a given environment U(x) the probability distribution for the position of a 
particle P(x,t\xo,0) satisfies the Fokker-Planck equation (in the continuum): 

d t P(x, t\x , 0) = d x (d x + U'(x))P(x, t\x , 0)) = -H FP P(x, t\x Q , 0) 

with the initial condition P(x, t — > 0|a;o, 0) — > 8{x — Xq). As is well known, setting 

G(X, t\x , 0) = e (U(x)-U{x ))/2 p ^ x ^ ^ Q) (227) 

one obtains the following imaginary time Schrodinger equation for the Green func- 
tion G(x, t\x , 0) 

d t g{x,t\x o ,0) = U% + \u"(x) - \u\xf\ G(x,t\x o ,0) = -H s g{x,t\x ,0) (228) 

with the initial condition g(x, t\xo, 0) — > 5(x — x ). This is the standard form for the 
Schrodinger operator H$ associated with a diffusion process. It can be factorized as 
H s = Q ] Q with 

Q = d x + U\x)/2 

and 

Qt = _ dx + U'(x)/2, 

and thus has a real positive spectrum E n . The Fokker Planck operator Hpp is non- 
hermitian but has the same real positive spectrum, with right and left eigenfunctions 
$^(x) and $n( x ) associated with E n . They are related to the eigenfunctions ip n (x) 
of the Schrodinger operator by &n( x ) = e~ u ^ x ^ 2 ^p n (x) and $^(a;) = e u ( x ^ 2 tp n (x). 

In the next two sections we use some of the results obtained previously for 
the Sinai diffusion process to obtain results for the Schrodinger and Fokker Planck 
operators. 



B. Averaged Green's function for the Schrodinger operator 

Interestingly, one can obtain the averaged Green's function of the Schrodinger 
operator ( [22 7|) from a slight variation of the previous calculation for the dynamics 
inside a well of Section |VI F| . The physical reason is that in Sinai's model the 
particle tends to jump to and occupy lower accessible wells, with weight e~ u ^ near 
the bottom. As a result one can show that, due to the the exponential factor in 
(227), the dominant contribution in the average over disorder of ( |227| ) comes from 



rare configurations in which the point x and the point are at about the same 
potential. The calculation is sketched in Appendix |E]. The result is 

WAt) = (£) ( 229 ) 

with r = hit with the scaling function G(X) given by (|1 74j) . 
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In the case of a small bias 5 > the result becomes: 



T5 \ 2 



Q{x > 0,0, t) = Q(x <0 1 0,t) = - I inh(r ^ j '*'G(|x|/r 2 ) (230) 

which is valid in the usual scaling regime T5 and x/T 2 fixed with T = lnt large. 
Note that this averaged Green's function Eq. fl230p , is closely related to the average 
Green's functions of a one-dimensional lattice fermion problem with random nearest 
neighbor hopping, t n = t + 5t n as computed recently by L. Balents and M.P.A. 
Fisher |j32"| . In particular, the inverse Laplace transform of Eq. ( |230| ), Q(x, xq,E) 
which is a function of the wave functions of 7i s at energy E, is equivalent to the 
Green's function of the Fermi problem at energy e = y~E with x — xq = n — no 
even, (corresponding to ipR + ipL m terms of the right and left moving fermions of 
ref. [|32|]). This is related, in the Sinai problem, to the dominance of the averaged 
Green's function by x and Xq both at bottoms of valleys which correspond to even 
sites. The random hopping, (— l) n 5t n corresponds to U'(x)/2 in the Sinai problem. 

By methods similar to those used in the present paper, one can obtain much more 
information about the statistical properties of eigenf unctions and Green's functions 



of both problems. These will be analyzed in ref. [48 . 



IX. CONCLUSIONS 

In this paper we have developed a powerful real space renormalization group 
(RSRG) method procedure for models of diffusion in one-dimensional random po- 
tentials which belong to the universality class of the Sinai model. This method is 
simple to implement, physically transparent and allows one to obtain exact results 
for universal quantities. 

The RSRG was first applied to recover, as a check of its validity, the single time 
diffusion front for the rescaled position x(t)/(T\nt) 2 obtained previously by Kesten 



and Golosov |12|,p7[ . In addition we obtained the diffusion front in presence of a 
small bias in the crossover region. 

The study of persistence properties, i.e., probabilities of return to the origin 
and their associated decay exponents, showed that in disordered systems distinc- 
tions must be made between recurrence properties of thermally averaged trajectories 
(x(t)} (exponent 9) and single particle trajectories (exponent 9). Nontrivial expo- 
nents (e.g. 9 = (3— \/5)/4) were obtained for thermally averaged trajectories, a novel 
and unexpected feature of the Sinai model. The distribution of number of returns k 
was found to be strongly peaked in the rescaled variable, g, at g = k/ ln(T In t) = 1/3 
but with multifractal tails characterized by an exponent 9(g). It was shown that 
single-run averages ~ f*x(t')dt' obey the same scaling for g < 1/3, but with devia- 
tions on the larger than typical (g > 1/3) side of the distribution due to rare events 
which were analyzed. We found that at each return to the origin, the thermally 
averaged trajectory loses some memory of the past. This allowed us to compute 
exactly the probability distribution of the complete sequence of return times. By 
contrast the successive jumps of (x(t)) exhibit persistent correlations which we have 
studied in detail. Much of the analysis was extended to the case of a small applied 
bias. 
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Single particle properties, such as return probabilities, distributions of first pas- 
sage times and of maximum displacement, were obtained by studying the RSRG in 
the presence of boundaries. The first meeting time distribution of two independent 
particles was also obtained. Extensions to large but finite size systems was studied, 
for reflecting or absorbing boundary conditions. The distributions of equilibration 
time and position, of the first passage times in presence of boundaries and the finite 
size diffusion front were all obtained. 

A second set of results concerned aging dynamics. The scaling form of the joint 
distribution of positions x(t') and x(t) at two times t' < t was obtained explicitly in 
Laplace transforms. This two-time diffusion front was found to possess an overall 
hit/ hit' scaling, with an interesting singularity at x(t) = x(t'). Explicit expressions 
of several moments and correlation functions were obtained. In the presence of a 
bias, our single time diffusion results (the distribution of x/t^ being related to a 
Levy distribution) and our two-time aging results (with a t/t' scaling) are consistent 
with known exact results and with the phenomenological description in terms of 
an effective directed model with an algebraic distribution of waiting times. But in 
addition we have obtained the full crossover between the symmetric and biased aging 
scaling forms. Our aging results are also consistent with the numerical simulations 



and qualitative arguments of [16 



We have also obtained several quantities which are controlled by rare events such 
as renormalized valleys with degenerate minima or degenerate barriers. These can 
be studied systematically as subdominant contributions in the RSRG. From them, 
we computed the fluctuations in the thermal width of the single time diffusion front 
(i.e., moments such as (x 2 (t)) — (x(t)) 2 ~ T(Tlnt) 3 ), the two-time diffusion front 
in the quasi-equilibrium regime (for t — if ~ t' a , a < 1). 

This work exhibits the relationships which exist between the Sinai model and 
problems such as quantum spin chains with disorder: both can be treated via very 
similar RSRG methods. Although observables of physical interest are often different 
in each of these models, some interesting connections, have appeared — e.g., between 
persistence properties of Sinai model and magnetization in the random transverse 
field Ising model. The RSRG methods enable one to consider this class of models 
in a unified way. Since the method allows one to check its own range of validity, it 
may shed light on different universality classes. The averaged imaginary time Green 
function of a related random Schrodinger problem was found as a side benefit. 

In conclusion, the model studied here provides an all too rare explicit example of 
a zero temperature glassy fixed point where detailed non-equilibrium quantities can 
be obtained. Qualitatively similar behavior should be expected in systems where, as 
in Sinai's model, the dynamics consists of jumps over large barriers between partially 
equilibrated configurations. The detailed understanding of physics in the simple one- 
dimensional case studied here perhaps encourages hope that new methods — exact 
or approximate — based on similar ideas can be developed for more complex glassy 
systems. As a start, we have already applied the methods introduced here to more 
complex one- dimensional models, in particular the non-equilibrium dynamics of the 
classical random field Ising model ||30|| , as well as reaction diffusion models with 



disorder Furthermore, recent work on random quantum Ising models in two 

and three dimension |^9| suggest that in at least some systems the type of behavior 
found here are not particular to one dimension. 
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APPENDIX A: AUXILIARY VARIABLE RG RULE, SYMMETRIC CASE 

In this Appendix we study the general RG rule m' = ami + brri2 + cm^ upon dec- 
imation of link 2 (see Fig. [I]). We introduce the rescaled variable fi = m/Y^ where 
ijj is an unkown exponent and look for the fixed point joint probability distribution 
P(rj, //) that is a solution of 

= ((1 + rj)d v + 1 + V (/i<9 M + 1)) P( V , AO 

poo POO POO 

+ djJLx I dfi 2 dfi 3 P(0,fi 2 )P{-,fii) * v P(.,/x 3 )5 (ju - {am + b/i 2 + c/i 3 )) (Al) 



We have of course P{rf) = J °° dfiP(r],n) = e v . The equation for the first moment 
C ( 1 l) = pfe) io°° d W p (v, M) reads 



= (1 + v)d v C(ri) -(ti + ^)C( V ) + (a + c) / d^C^) + bC(0)r) (A2) 

It is useful to differentiate this equation with respect to rj to obtain 

= (1 + v)d 2 v C(ri) + (l-v~ ip)d v C(rj) + (a + c - 1)C( V ) + bC(0) (A3) 

and to keep the boundary condition C"(0) = C(0)ip at r] = 0. 

For a + c—1 ^ (the case for all the physical quantities discussed in this paper), 
it is convenient to set y = 1 + 77 and T{y) = C{rj) + a+ ^_ 1 C(0) so that that T(y) 
satisfies the confluent hypergeometric equation 

= yd 2 y T +{B- y)d y T - AT(y) (A4) 

where B = 2 — ip and A = 1 — a — c, together with the boundary condition T'(y = 
1) = ip-rzuT(y = 1). Since we are looking for a well-behaved (i.e., not exponentially 
growing) solution at r] = 00, we find that T(y) has to be proportional to the confluent 
hypergeometric function U(A, B, y). To satisfy the boundary condition at y — 1, we 
obtain, using the functional relation U'(A, B, 1) = U (A, B, 1) — U (A, 5 + 1,1), that 
the exponent ip has to satisfy the equation 

= ( 1 - i\) a + c ~ l \ u(l-a-c,2-il),l)-U(l-a-c,3-il>, 1) (A5) 
\ a + o + c—lj 

Note that in the particular case where A = 1 — a — c = — 1, the function U(—l,B,y) 
reduces to the linear function — B + y, and the equation for ip is simply quadratic 
ijj(ip - 1) = 1 + b yielding ip = (1 +_y/5 + 4b)/2 as in 

In the text we use the ratio m/lr which decays as ui/Iy ~ with $ = 2 — ip. 
Both exponents ip and $ depend explicitly on the coefficients a, b, c, of the RG rule. 
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APPENDIX B: AUXILIARY VARIABLE RG RULE, BIASED CASE 

We consider the auxiliary variables m ± that evolve with the RG rules m + = 
a + m± + b + rri2 + c + mf upon decimation of an ascending link 2 and vrC = a~m\ + 
b~m\ + c~m^ upon decimation of a descending link 2. We introduce the joint 
probability distributions P^((,m) that evolve as 

POO 

(9 r -9 c )P r ± (C,m) = P r ± (C,m) / dm 2 (P± (0 , m 2 ) - Pjf(0,m 2 )) 

Jo 

/*oo /*oo /*oo 

+ / dm\ I dm 2 / dm^P^ (0, 7712) 
Jo Jo Jo 

P^(., mi) *^ P^(-, "^3)5 {jn — {ami + b ± m 2 + c m 3 )) . (Bl) 

We have Pf^(C) = Jq 00 (imP^ (£, m) = u^e - ^ . The equation for the first moments 
Cr (C) = Jo°° rfm mP v (C, m ) is 



(ar-a c )c±(c) = 4% 



z (&±^(0)-C±(C)) + (a ± + c ± ) / dC'C^C) 



(B2) 



with u+ % = = 1/I r 

We study the simpler particular case when a + + c + = 2 , a~ + c~ = 2 , b + = 
b~ = b. Then the solutions C^(() = C^(() = Ct(C) are simply linear in ( : 
Cr(C) = ^4r + (Br and the coefficients satisfy 

B r = d r A r (B3) 

d 2 v A T = (1 + b)u^A T = (1 + b) —^— A v (B4) 

For 5 = we have already seen in Appendix [A| that the auxiliary variable m 
grows as T^ with ^(6) = (1 + v / 5 + 4&)/2. Indeed for 5 = 0, A r oc is a 

solution of (|B4]). For 5 > following [^TJ we thus look for a solution of the scaling 
form: 

A T = r^ (f> M (fe) (7 = ST) (B5) 
where A^^) satisfies the equation: 

d^(j) = ^^ ( %) (B6) 

with the boundary condition A^ip/) oc 7^(6) as 7 — > 0. Introducing the new variable 
y = coth7, we obtain the differential equation for the Legendre functions: 

d^Mb) j Mb) 

The solution for ^^(7) with the above boundary condition is: 
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A^ b \ 1 )=KQ m . 1 (coth 1 ) 

= K> tanh( 7 )^F (M±± m + 1, tanh( 7 ) 2 ) (B8) 

with fT' = KypH2~^^ i r ^h)+i/2) w h ere is a non-universal constant. The 
asymptotic behaviors are A^(7) ~ K'^^ at small 7, and ^^(7) = 
2if ' p7?^v| i)/2]r[^(b)/2] 7 a ^ l ar S e 7- We can now compute the mean values of the 
variables m 

poo poo poo Q A 

(m ± )=/ dC / dmmP±(C,m) d(P±(()C±(() = A r + ^ (B9) 
Jo Jo Jo M r 

yielding fl92|) in the text. 

APPENDIX C: CORRELATION OF TIMES AND DIRECTIONS OF 

SUCCESSIVE JUMPS 

1. Conditional probabilities of times of jumps forward and backward 

In this Appendix we compute the conditional probabilities pp r , to make a jump 
forward at F (respectively a jump backward) given that the last jump occurred at 
T'. We define Dr t r'{F) as the probability to be on a descending bond of barrier F 
given that the last jump of the effective dynamics occurred at V (this jump was 
necessarily in the (+) direction since the walker is on a descending bond). The 
initial condition is thus given by: 

POO POO POO POO POO 

D rr (F) = K dF 1 j dh dl 2 dF 3 dhih + h) 

P V {F X , h)Pr(T, k)P r (F 3 , l 3 )5 (F - (F, + F 3 - F)) (CI) 

Indeed, the bond must be a new bond made, at T', out of three bonds, and the origin 
of the random walk must have been on either the first or the second bond in order 
to satisfy the condition that the last jump occurred at scale V . The normalization 
constant K has to be choosen to ensure that f T Dt',v(F) = 1. Introducing the 
rescaled variables 77 = and a = ^7, we obtain that D a (rj) evolves according to 

(ad a - (1 + T])d n - 1) D a {rj) = -2D a (r)) + f dr/D a (r/) e -^ (C2) 

Jo 

with the initial condition at a = 1 given from ( |Cl~D by D a= i{r]) = {^ + e~ v . 
The solution reads 

D a ( V ) = A a e~^ + (B a + C aV )e- a ^ (C3) 

I - ' ■■ - ( C4 ) 



.4 



2a 



2 



1 / \ 2 
If a 



B « = -2{^l) ^ (C5) 

C a = -U -At! (C6) 



2 \a 



65 



The probability to make a forward jump at T (i.e., in the (+) direction) given that 
the last jump occured at r" (and by convention was in the (+) direction) is 

/OO 1 /»oo 

dFD T>T ,(F) = - J d V D a ( V ) (C7) 

since the probability that the neighboring bond is decimated at T is -Pr(r) = 1/r. 
On the other hand the probability to make a backward jump at T (i.e., in the (-) 
direction) given that the last jump occurred at T' is 

f$v = DrAr) = f D °M ( C8 ) 

which is the probability to decimate the bond we are interested in. The total prob- 
ability to jump at T in any direction given that the last jump occurred at T' must 
satisfy pr,r' = Prr' + Prv = dFD rr <(F). These expressions, after substi- 

tuting the above solution ( |C6D yield the formulae (^) given in the text. 



2. Correlations in the sequence of times of successive forward and backward 

jumps 

A full calculation of all terms is quite involved and goes beyond the present 
paper. Here we indicate only the result for the two first elementary building blocks 
for the many jump correlations. The first one is 

p^T^d^ = p^{a l )da l = ^ (2 - (1 + a l )e-^) , (C9) 

which is a scaling function of ai = Fx/To- Intermediate calculations also yield the 
probability that the second jump occurs at Ti and is a forward jump given that the 
first one occurs at T and is backward. 

^(iMro^r! = pi\a 1 )da l = (2 - e"^- 1 )) . (CIO) 
The second elementary building block is given by: 
Wfr. „r.ur^. _ ^ ~. ^- - _ da,da 2 (2 - e~^) 



p{ (T2,Ti\T )dT 2 dTi = p{ (a 1 ,a 2 )da 1 da 2 



a\a\{2 - e-^i" 1 )) 

-(a2+ai-2) p -(a2ai-l)~ 



( -(o 2 +ai-2) p -(a2Qi-l)\ 
4 - e"^" 1 ) - 2(a 2 + l)^" 1 ) - (a 2 + 1)- — + {a 2 a x + 1)- — (Cll) 
Qtl — 1 «l — 1 / 

which is a scaling function of a± = Ti/Tq and a 2 = T 2 /Ti. 



APPENDIX D: DISTRIBUTION OF SEQUENCES OF RETURNS TO 

THE ORIGIN: BIASED CASE 

To compute the conditional probabilities p ± (T, V) of returns to the origin defined 
in the text (Section |IVF|) we consider the probability D rr ,(£) that a bond has barrier 
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z at T and has not changed orientation since the scale T' where it was created. Its 
RG equation is 

(dr - C ) L»± r ,(C) = 2Pjf (0)P>± r ,(.) H P±(.) - 2P r T(0)L»± r ,(C) (Dl) 

with the initial condition Dp, r ,(£) = («r/) 2 Ce -1 *" 1 ". Since P^(C) = «pe"^r, the 
solution has the following form 

L>± (C) = (4^ + C5r,r0 K)" e"^ (D2) 



where the coefficients are 



A T> r> = 5~ l (g^i(coth 7 )P^ 1 (coth 7 / ) - g _ 1 (coth 7 / )P^i(coth 7 )) (D3) 

B t ,t> = —^2- (P;_ 1 (coth 7 )g^ 1 (coth 7 / ) - P _ 1 (coth 7 ')g;_ 1 (coth 7 )) (D4) 
smh 7 

where 7 = 5T, 7' = 8V, (ft = (l + \/5)/2 and Q v (y) and P,(2/) are associated Legendre 
functions: they are two linearly independent solutions of the equation (|B7|) (with 
(1 + 6) -> v{l + v)). 

The probability for a (±) link to be decimated at V given that its last decimation 
occurred at V is therefore 

p ± (r, r') = -dr / dCDfr, (0 = ^4A rr (D5) 
Jo u v 



This leads to the equation (p. 01] ) given in the text in terms of the the reduced 



variables y = coth 7 and y' = coth 7 '. 

APPENDIX E: DYNAMICS WITHIN A WELL 

1. Probability that a bond has degenerate minima 

Let us introduce the probability Sr(C> ^ x > w ) that: a given point (denoted x in 
Fig. [J) belongs at T to a bond of barrier F = V + £, of length Z, is at a distance 
x from the min of the bond and is at a potential w above the potential of the 
minimum of the bond. One has that by definition < x < I and its normalization 
with respect to x and w is J^dx J Q r+ ^ dwSr((, I, x, w) = lPr((,l)/ J, Pr(0> which is 
the probability that a given point belongs to a bond with P, I. The RG equation for 
Sr(C> 1} x i w ) reads: 

(dr - d ( )S r (C, I, x, w) = -2P r (0)S r (C, I, x, w) + P r (0, .) P r (., .) * cl S T (., .,x,w) 
Sr((x,h,x 1 ,w 1 )Pr(0 1 l 2 )Pr(C3,h) 



' Gi,h,xi,wi,l2,G3,h 

S(C - (Ci + Ca))S(l - (h + h + h))5(x - (xi + Z 2 + h))5(w - (w, + Ca)) 



+ / Pr(Ci,h)Sr(0,l 2 ,x 2 ,w 2 )P r (C3,h) 

S(C - (Ci + Cs))5(l - (h + h + h))5(x - (h +l 2 - x 2 ))5{w - (w 2 + Ci)) (El) 
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(a) (b) (c) 




FIG. 9. Different terms which contribute to the RG equation for S. 

where each term is described in Fig. (|9|) We notice that the evolution equation for 
Sr(C) l,x,w — 0) decouples and leads to the form: 

S r (C, l,x,w = 0) = -^(5(x)P T ((, I) + R r (C, I, x)) (E2) 

where the delta-function part represents the probability that the point Xq happens 
to be exactly at the bottom of the renormalized bond in which case w = by 
definition. The function Ry(C,1,x) is the probability that a renormalized bond at 
scale T has (£, I) and a distinct degenerate minimum at a finite distance x. For small 
x (x = 0(1)) this function is non- universal. We compute this function in the scaling 
regime x ~ T 2 where it is universal and of order 1/r. We use the rescaled variables 
rj — C/r, A = l/Y 2 and X = x/T 2 such that P r (C I, x) = r" 6 P(r/, A, X) and obtain 
the following fixed point RG equation for R: 

= Td r R r (r], A, X) = ((1 + V )d v + 2Xd x + 2Xd x + 6)R(r), A, X) 

+P(0, .) * A P(., .) *„, A P(, ., X) + 2P{ V , A - X)(P{0, .) * x=x P(0, .)) (E3) 

This equation was obtained by substituting the decomposition (|E2|) into (|E1|) in 
the spirit of an expansion in powers of 1/r, where the zeroth-order equation is 
satisfied by the delta-function part. The order 0(1/T) equation yields the equation 
for Rr((,l,x) where the delta-function part acts now as a source in the last two 
terms of ( |E1|) leading to the last term in (|E3|) . This term describes the probability 
that between T and r + dT a new bond with one degenerate minimum (a distance 
of order T 2 away from the lowest edge) is created via the decimation of a bond 
whose neighbor also has ( ~ — cases (b) and (c) in Fig. [| with w = w\ = 0) — the 
probability of this is of order 1/r. 

Before proceeding further, we notice that it is easy to also keep track of the barrier 
r between the two degenerate minima. We define Rr((, I, x, Tq) as the probability 
that a renormalized bond at scale T has ((, I) and a distinct degenerate minimum 
at a finite distance x separated from the minimum by a barrier T . It takes the 
scaling form Ry((,1,x,T q ) = r _7 i?p(?7, A, X, u) with u = T /T. The normalization 
is Rr(r],X, X) = J" 1 duRr(r}, A, X, u). The scaling form satisfies the fixed point RG 
equation: 

= Td r R(v, A, X, u) = ((1 + r])d v + 2Xd x + 2Xd x + ud u + 7)R(r], A, X, u) 

+P(0, .) * A P(., .) *„ >A R(., .,X,u) + 2P( V , A - X) (P(0, .) * x=x P(0, .))5(u - 1) (E4) 

where the last term corresponds to barriers T = T created upon decimation. 

Remarkably, one can find the complete solution of this equation in a factorized 
form: 
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R(r), A, X, u) = P(r], A — X)r(X, u) (E5) 
where r(X, u) satisfies 

= (2Xd x + ud u + 4)r(X, u) + 2(P(0, .) * x=x P(0, .))<*(« - 1) (E6) 
whose solution is: 

r(X, u) = 6(1 - u)u- 4 f(X/u 2 ) (E7) 

f(X) = 2(P(0,.)*A=xP(0,.)) (E8) 

which, using (|TT|), yields formula (|170|) given in the text. 

The biased case can be studied similarly. The corresponding quantities (as 
usual ± designate descending and ascending bonds, respectively) also satisfy 
S^((,l, x,w = 0) = j-(5(x)P^ + K^((,l, x)) and one finds that the RG equation for 

Rt(CJ,x,r ) is 

(dr - d c )B4(C, I, x, T ) = (P r ± (0) - Pjf (0))P±(P, /, x, T ) 

+P r ^(0, .) P±(, .) * cl P±(., ., x, T ) + 2P r ± (C, I - x)P?(0, .) * x P r ± (0, .)5(T - T ) (E9) 

The solution again factorizes into 

P±(C, /, x, To) = Pf((, I - x)rr(x, T ) (E10) 

where r r (x, r ) = 2#(r — r )P r t ) (0, .) * x Pr~ o (0, .) does not depend on the direction 
of the bias. Its Laplace transform is simply 

f°° r> + 5 2 

/ dxe-^r r (x,T ) = 9(T - T )2U^(p)U^ (p) = 29(T - T )—j^-^==- (Ell) 
Jo smh (r VP + ) 

so that finally 

r r (x, To) = 9(T - T )±r ( e-* (E12) 



F 4 \T 2 
1 o V 1 o 



where f is the function for the symmetric case introduced in (|170|) . 



2. Relationship to the associated Schrodinger operator Green's function 

The disorder averaged Green function defined in Section |V111 B| is exactly related 
to the probability Sr((,l,x,w) introduced above. In the symmetric case, one can 
restrict to x > 0, and one has: 

i roc K+T roo 

G(x,0,t) = - d( dwe~ w/2 dlS r (C,l,x,w) (E13) 
2 Jo Jo Jx 

The factor ~ is simply the probability to be on a descending bond [x > 0). This 
can be expressed using the rescaled variables ( = Trj , w = Tu , I = T 2 X , x = T 2 X 
and simplified using that for large T, we may replace e~ Vu l 2 by ^8(u). Using 
one obtains: 
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2 poo roc 

G(x,0,t) = — drj dX R(r},\X). (E14) 
1 Jo J x 



Using (p5|) one finds the result of the text fl229Q . 



In the biased case we obtain an expression for the averaged Green function in 
terms of the functions Sp((, Z, x, 0) 

poo poo 

G{x > 0,0, t) = 2 d( dlS£((,l,\x\,Q) (E15) 

J0 J\x\ 

rOO POO 

G{x < 0,0, t) = 2 d( dlSr(F,l,\x\,0) (E16) 



Using ( FTUD we finally get £0 > 0,0, t) = £(x < 0,0, t) = G(\x\,t) with 

G(\x\,t) = ^- £ dV r r (\x\,V Q ) = lit? fM) e -l^ 2 (E l 7 ) 

where the function G has been introduced in ( |174j) , leading to formula ( |23C| ) in the 
text. 

APPENDIX F: SOLUTION OF THE TWO TIME RG EQUATIONS 



In this Appendix we solve explicitly the RG equations ( 147 ) for the quantities 



fip r ,(£, Xl, xr; x' l , x' r ), e = ±1 e' = ±1 from which we can obtain the two-time 

diffusion front Prob(a;t, x't'\00). We consider the general biased case and discuss the 
particular limit of the symmetric case. 

We first introduce the Laplace transforms: 

= [ e-^ v «e-^e-^ fl -^ ) e- v(si -^ ) n^ r ,(C,XL > Xfl;x£ > x / Ji ) (Fl) 

Jx' >0,a;' D >0,3;r,>0,a:R>0 



Since we consider large T' we can use the fixed point solution 



(F2) 



poo 

P r ± (C, / u)= / rf/e-^P r ± (C,/) = f/ r ± (/i)e-^r(M) 
Jo 

The RG equations (|147|) can be then written in Laplace variables as: 
(dr - d ( ) Cl^, (C, fi, v- //, v') = -2£^f(0)fi±£ (C, fi, u- //, v 1 ) 

POO 

Jo 

poo 

+U+(v)U r (v) / dC^^Mfl&iCxiWpW) 
Jo 

poo 

+n^,(0,fi,is^',is')U±(fi)U±(v) / d&e-U-^Me-^M (F3) 

Jo 
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together with the initial conditions at T = V given in (|148|) which become: 

/>oo /"OO 

fif^r, (C, li, v; //, i/) =5 ee , / dx' L / dx' R e-^ x '^e- v ' x '^u( l {^x' w x , L ) 

Jo Jo 

= S^-^—f (U^v')e-^ - U^')e-^>'A . (F4) 
'r' ^ 

We look for solutions of the form 
Cl^ r ((, /x, i/; //, v>) = Af[ T , (ji, v- //, v') e-WM + „ ; ^ ,/) e -d*M 

+C^v (/i, I/; //, l/) e^rV) + ( At , z/; //, */) e"^r' ("'). ( F5 ) 

It is useful to introduce the functions 
ef[ v , (fu, v\ fi', v') = £}f e ' r , (C = 0, fi, v\ //, v') 



A v,v (p> v\ /A v ') + B r,r> (p> v \ tf> v ') + c r,v (j*, v \ /A v ') + D r,r> (p, "5 li', v) 



ofjv (ll, v, ll', v) = / d£Qrjv ((, li, v\ //, v) 
Jo 



_ Af'r (p, v\ li', v') + ggp {fr v _\ fA v ') + Cf'r/ (/x, v\ ll', v') + gf^v (//, ^ //, z/) 
wf(/x) U \W) Ur'd^') uf(u') 

For each initial condition indexed by one of e' = ±1 it is most convenient to work 
with the eight functions 8 ee , a ee , C ee and D ee with e = ±1. 

We first consider the equations for C ee and D ee . These equations are homoge- 
nous and thus easier to solve. The equations for C ee read 

drCf r , (li, v; li', v') = ( {0)6 e , +1 - 2*7+(0)5 e ,_ 1 - v&(jjl') + Ur ( ^ )f/r+(/i) 



Ut(li) -u^(li') 



«f(V) - «f/(/i') 
with initial conditions at Y = V 1 indexed by e': 



C^v (ll, v; J, v') = -«y * 3^1. (F8) 

tp/ Ll V 



Similarly the equations for D ± read 



drDf r , (li, v; ll', v') = ( -2£7f (0)<J 6i+ i - 2[/+(0)5 e _i - ufV 



Uf(n)U£(ji) 
uf(/x) - wf',(z/) 



^J£p^Mi D 0^ u . M {F9) 



with initial conditions at T = V are indexed by e' 



D^ r (/i, i/; ^, iO = (F10) 
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To find the solution one notices that each matrix element C te ' and D ee ' satisfies 
its own differential equation. Thus, since the initial condition is diagonal in ee', the 
solution is also diagonal. It is found to be: 

Cr'r> fa v\ n\ v) = 

u T (Q) 2 u^fa) (4fa) -^'M)KV) - 4^)) c - ( r-r'w;, w 
i T , uf(oy hi' - 1/) - i4(ji))(i4,(ji') - i4(v)) 

D r'r> fa v\ fi', v 1 ) = 

5* tff (0) 2 ^(V) KV) - ^,(/i))KV) - c _ ir _ r , w > Al/l) 

l T , Uf(0) 2 (// - v') (u<t,(v>) - i4(ji))(i4(i/) - u^{v)) 1 ' 

For each e' the four remaining functions 9 ±,e ' and a ±,e ' satisfy the following 
system of four differential equations: 

d r d +t ' = -(2u r (0) + u+(fi) + u+(u))9 +€ ' + u^(fi)u+(u)a +e ' + p+ r 5 e ', + (F12) 
<9 r ^ e ' = -(2u^{0) + u r {fi) + u r {u))d- €> + u r (n)u r (u)a- e ' + p^ r ,5 ef ^ (F13) 

= f-2„ r - ( o) + + m^m ^ _ ^ ^woy (F14) 

= f_2 Mf (o) + + ^MgM) _ r , + (F15) 

We note that the system for e' — +1 and the system for e' = — 1 are identical 
except for the inhomogeneous part, that we have defined as: 

Pr,r' = — TT^farfa ~ «r'(^))(«rV) - u r>fa)) c r$> 

U r ,{fi ) 

-^—AuHn) - uUu')){4{u) - 4,(v'))D&. (F16) 
To exhibit explicitly the F dependence, it is useful to rewrite 

p+ r , (/,, v- //, v>) = 1 (f+fa v- M ')e-( r - r ')v(M') - /+(,,, „ ; ^ )e -Cr-r0v<'>) (F17) 
It\ji'-V) 



p rr fa i/; //, i/) = ttt M r fa "I //K (r - r »' } - fcfa „ ; ^e^r')^')) (F18) 
where we have introduced the two T-independent functions 

/r'< fa v\ /0 = ^T7V(4''(^) - t£(A*))(t£(A0 - (F19) 

The first step is to obtain the solutions of the above two (identical) homogeneous 
systems ( F15| ). Remarkably, these can be constructed from the functions U^fa) and 
v$fa) using the differential equations (ED). We find the four independent solutions 
of the homogeneous system to be: 
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n r U^(fi), d 1 = -n T U r (u) 



-n T 



(F20) 



nrCp (/x)wp(i'), 9 2 = n r U r (v)u T (//), <r 2 



n r - 



x «(^-f/ r »[/ r >)) ; 



(7 9 



n r 



6+ = n r U T 



3 = -n r U r (fx), a 3 



(ur(Li) 2 -U+(Li)Ur(Li)) 
WW 



(J, 



u T [y) 



(F21) 
(F22) 



Q~\ = n r Uf (u)u^(fi), 6 4 = n r U r (li)u t (u), o\ 

e£-(aO 



Ur 



nr 



Up (^)Wp (v) 



u v {ll)u v (u) 



X( W +(/i) 2 - U+((J,)U T (//)), (74 

where n r = l// r = sinh f (r ^ and with <9 r n r 
eful to consid 

/ 



-U+(u)Uf(u)) 



(F23) 



Mr + u v )nr. 



sinh a (r\S) VV1U11 ~ _ V"T -r "T, 

It will be useful to consider the matrix formed by these solutions 



Nr 



Qz 6 



\ 



(To 



'1 °2 "3 u 4 

From the usual properties of systems of linear equations the Wronskian Wr = 
det[iVr] satisfies the simple equation, drWr = Tr[Mr]Wr where M is the matrix 
formed by the coefficients of the homogeneous part of the linear system. One can 
easily integrate this equation, or one can compute directly the determinant, and use 
the definitions ( F23|) to simplify the result (after a tedious calculation). Both give 
the same, remarkably simple, result: 

^U+(li)U^(li)U+(p)U^(u) 



a 



0~-i 



\ 



J 



Wr 



-n r (/i 



(F24) 



Since this is not zero, this shows that the four solutions given above are linearly 
independent. Thus we are now in a position to write the solutions of the two linear 
differential systems ( [F15| ) with the inhomogeneous terms. It is found, as usual, as 
a linear combination of the four independent solutions (|F23|) of the homogeneous 
system: 

/ &trM \ ( ^ 

af r (Li, v) 

App/(/i, v, //, v') are the coefficients of the linear combinations. Using 
the standard method one finds the following equations for the coefficients: 

/ d r X 1+ \ ( vtr* \ 



( 



\ 



+c' 



/ 



i=l 



Nr 



w 



\ 



/ 



(F25) 



where Ap 



iVr 



2+ 



<9 r A 3 + 









(F26) 
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iVr 



d r X 2 - 
\ <9 r A 4 - 










\ 






Pt,t' 















V 





/ 



The initial condition for the Ap 6 ' at T — V are fixed by the initial conditions 

4 

6r',r> (/■*> v\ v ') = E ^' ^ °^ T ' J ^ ^ = ^'f (C = 0, fi, v\ n', v') 
i=i 

= 6ee'J--^— f fcV) " ^V)) (F27) 
<p/ /i — V \ / 

^ poo 

0r''r' (/i, z/; //, i/') = A^'p 0, v\ /A ^) <(r', /i, i/) = / rfC^F''r' (C> M> "I /A ^) 

i=l - 70 

-5 1 1 I^V) U^')\ (F28) 

Finally, we find, for e' = +1, the solution of ( |F15|) with the above initial conditions: 
4 r „ r 



i=i \- Jr ' 
4 r r 

°r}> = E a tr [ dT(Nr\ lP + r/ + (N^) i , l e++, + (N~\ 3 a^ r , 
i=i 



(F29) 
(F30) 



and, for e' = — 1: 
l 



i=i L^r' 
4 r ,r 

i=i \- Jr ' 



(F31) 
(F32) 



The next step is to evaluate iVjT , the inverse of the matrix Nr- Remarkablythe 
inverse admits a simple explicit form in terms of the functions and U^r, which 
can be found after some tedious calculations using the form (|i9|). It reads: 



n T {v - ji) 



A* 


V 




tip (u)itp (i/)u r (i/) 




Uf(u) 




t/p» 




u r( u ) 




Up (/i)«p (l/) 






[7+00 




V 




«P (u)up~ (^)up (^) 


Up (m)«p (a*) 


t/p» 




c/p+H 




u+(u) 




M+(([i)w+(y) 


Up 


U+(u) 


Up-( M ) 


£#(«/) 





(F33) 



We have thus obtained the quantities of interest, namely the crf e r /(/z, v\ fi'u') = 
C f2 r e ' r ,(/i, z/; /iV). The two-time probability can then be obtained from the a as 
follows. 
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The Laplace transforms in the different sectors (x > 0, x' > 0), (x < 0, x' < 
0),(x < and x' > 0) and (x > and x' < 0) are respectively: 

poo poo 

P++(p,p') = dx dx'e- p ' x 'e- p ( x - x,) Prob(xt,x't'\00) 
Jo Jo 

= a++, (p = p, v = 0; p! = p', v' = 0) (F34) 

/0 pO 
dx / dx'e p ' x 'e p{x - x,) Pioh(xt,x't'\00) 
■oo J — oo 

= crp p/(/x = 0, v = p; // = 0, v' = p') (F35) 
/ dx / e -p ' x 'e p(x - x ' ) Prob(ar*,x't / |00) 

oo JO 

= o-p p ; (/i = 0, f = p; fj! = p + p', v = p) (F36) 

/0 poo 
dx' / rfxe p ' x 'e- p(x - z,) P7ob(^~?Fp) 
■oo </0 

= dp p, (// = p, v = 0; fi' = p, v' = p + p) (F37) 

From these one can compute the distribution Q(y,t,t') 
= f_°° dx'Prob((x' + y), t, x', t'\0, 0) of the relative deplacement y = x(t) — x{t'). 
Its Laplace transform in the sectors y > and y < are respectively: 

POO 

Q$Ap)= / dye-™Q{y,t,t') 
Jo 

poo poo 

= dx' dxe- p{x - x,) Piob(xt,x't'\00) 

Jo Jx' 

/0 poo 
dx \ dxe- p{x - x,) Piob{xt,x't'\m) 
■oo Jo 

= o-p p,(/i = p, v = 0; // = 0, i/ = 0) + o^p/t^ = 0, v — p; fj! — p, u' — p) (F38) 
Qt,Ap)= f dye p yQ(y,t,t') 

J —oo 

/0 p0 
dx' / dxe p{x " x,) Prob(xt,x't' |00) 
oo J — oo 

/0 ^oo 
cix / dx'e p{x ~ x,) PToh(xt,x't'\00) 
-oo Jo 

= CTpp,(/i = 0, = p; ji = 0, i/' = 0) + erf p,(/i = p, v = 0; // = p, v' = p) (F39) 

In order to compute these distributions, we have explicitly evaluated the sums 
and integrals in ( |F30| , |F32| ) . As an example let us examine: 

4 f F 

<7jtJ,(//,z/;f/,i/) = YV;r(/^)[/ dTiN^i^^pf 

i=i Jr ' 

+(iVf/) i)1 (/i, ^+ r ,(//, i/) + (A^M/i, ^^T'O"'* i/)] (F40) 

It turns out that all integrals that appear in this expression are simple exponentials. 
This remarkable property remains true for all other elements and is the reason why 
the calculation, though tedious, can be carried out explicitly for this problem. 
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Explicit results for the symmetric case 5 = 



Here we give the explicit expression for the Laplace transform of Probst, x't'\00) 
in the different sectors defined in ( |F37| ). For the sector where the product xx' > 0, 
we find: 

~ ++ / s 2 , r 2 \ ~ / s 2 , r 2 

"r= a r',r' I P = Y^ ,p = F 2 / = r=ar '. r ' \ p = p2' p = p2 



coths 
a 2 r 2 s 



1 - 



1 



cosh r) a 2 r 2 s cosh(as) 
cosh s coth s + ( 1 — — ) sinh s 



coshr 



sinhr 

(r coth r — 1 ) (r coth r — s coth s 
cosh r(s 2 — r 2 coth 2 r) 



■{— s cosh s — r cothr sinh s 



+e 



-((X—l)r coth r , 



s cosh(as) + r cothr sinh(as))}] 



— 1 + s coth s ( 1 



tanh(as) 

a 2 r 2 s |_ srnnr V coshr 

(r coth r — 1) (r coth r — s coth s) tanh r 



r cosh r 

For the sector where xx' < 0, we find: 



-y g— (a— X)r cothr^ 



(F41) 



■P 



r 2 — s 2 



r=ar',r' l /' 
1 

a 2 (s 2 — r 2 )s \ \ coshr 
1 



Y'2 ' ^ Y' 2 
( coshs 



dx'e r r ,: 



+oo 
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I ) coth 6 
1 



a 2 (s 2 — r 2 )s cosh(as) 



sinh s 



sinh s 
sinhr 
cosh s 
coshr 



(r cothr — l)(r cothr — s coths) 
cosh r(s 2 — r 2 coth 2 r) 



{—s cosh s — r cothr sinh s 



+ e 



-(a— l)r coth r , 



s cosh(as) + r cothr sinh(as)) 



tanh(as) 



a 2 (s 2 — r 2 )s 



cosh s 



+ 1 - 



cosh s 



sinh r sinh s 
(r coth r — 1) (r coth r — s coth s 
cosh r(s 2 — r 2 coth 2 r) 

- e- (Q - 1)rcothr (rcothrcosh(as) + ssinh(as))} 



coshr 

{r coth r cosh s + s sinh s 



(F42) 



and the same expression for P h . 

We have performed a similar calculation in the biased case, but the corresponding 
full expression for the Laplace transform is too lengthy to give here. Some particular 
limits are discussed in the text. We give here the explicit expression of the two time 
correlator: 



(x(t)x(t>)) 



1 



32<5 4 sinh^ i sinh 2 7 



A(7,7 / ) + 5( 7 ,7 , )e- (7 ~ 7 ' )coth(7 ' ) 



with: 
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A(7, Y) = cosh(27)(sinh(27) - 7)[sinh(4Y) + sinh(2Y) - 6Ycosh(2Y)] 



+ cosh(27) 



—7' sinh(47') + 



cosh (47') 



+ 4 7 ' 2 cosh(2Y) - 2Y sinh(2Y) 



1 



2 ' v ' ' ' " ' ' 2_ 

+ sinh(27)[— 7' cosh(47') + sinh(47') — 7' cosh(27') + sinh(27') + 27' — 67 /2 cotl^'] 
+7[-2Y(cosh(4Y) + 4cosh(2Y) - 2) + 3(sinh(4Y) + sinh(2Y)) - 12 7 ' 2 cothY] 

5' 

+27' sinh(67') — 2 cosh(67 ; ) 



6Y 2 + - ) cosh(4Y) + 9Ysinh(4Y) 



+ (6Y 2 - 2) cosh(27') + 4 7 ' sinh(2Y) - coth(Y) + 8Y + 



13 



(F43) 



5(7,7') = cosh(37 + Y)(sinh(2Y) - 2 7 ') 2 

+ cosh( 7 + Y) (sinh(2Y) - 2 7 ') [-87' cosh(27') + 5 sinh(2Y) + 4(7 - i) 
-477' cothY - 27'] + sinh(7 + 7 , )47 /2 (sinh(27 / ) - 27') 



7 



+ cosh7(sinh(27') — 27') 



sinh(37') + 47cosh(7') + sinh(7') — 47 



7' 



sinh 7' 



sinh(7') 2 

(F44) 



APPENDIX G: DISORDER AVERAGED PROBABILITY 
DISTRIBUTION FOR A FINITE SIZE SYSTEM 

In this Appendix we consider a finite size system < x < L using results of 
Sections |V11 A| and [II C 3 . We start with a RR system, i.e., reflecting boundaries on 
each end. 

We will denote by b = 2k + 2, k = 0, 1, 2, .. + 00, the number of renormalized 
bonds in the system. The disorder averaged distribution can be written as a sum: 



-00 2k+2 



P OL (x,t\x ,0) = 2^2fc+2( X >*|30,0) 



(Gl) 



fc=0 N=l 



where P^k^^-i ^l x o, 0) corresponds to the contribution of the case where the starting 
point xq is on the N th bond (see Figure |8[). One must distinguish between iV = 2n+l 
odd, when the particle starts on a descending bond and x > xq and iV = 2n + 2 
even when it starts on an ascending bond and x < xq. Thus, in addition to (|G1| ) 
above we will also be interested in the explicit decomposition: 



Pol(x, t\x , 0) = 6(x - x Q )P£ L {x, t\x , 0) + 6(x - x)P 0L (x, t\x , 0) 

fc=+oo k 

6(x -x )Pj~ L (x,t\x ,0) = Y ^2 P 2k+2' L ( x ^\ x o^) 

k=0 n=0 



One has (see Figure §), for n = 0, ..k: 

2n+l 

D 2k+2' { X ,t\ X 0,V) 



(G2) 
(G3) 



2n 



X 



(G4) 



,i=i 



2k+2 
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where (..)2k+2 denotes the average over the 2k + 2 bond sector of the finite size 
measure RR in ( 198 ). There is a similar formula for an even number of initial 
bonds. Throughout we will define Laplace transforms x — > p, x — > p , L — > q as 
follows: 

^2fc +2 (P,Po,g)= / di/ <*r/ dx e-^ + ^ o+ ^P 2 ^ 2 (x,t|x ,0). (G5) 
Jo Jo Jo 

One finds: 

-^2fe+2 (PiPo,q) — —E p+po+q (P p+po+q P p+po+q ) P p+po+q {P q+p — P p \ po+q )(P q P q ) E q 

^2k+2 (P,P0, (?) = —Ep + + (P + + Pp +po+q ) (P — P po+q )P q (P q P q ) E q 

Po 

the first formula being valid for 1 < n < k (and k > 1) and the second for < n < 
k — 1 (and k > 1). Finally for the two edge bonds one has: 

^(P.PO, ?) = - J?^)(P fl -P fl + ) fc JS7 (G8) 

Po 

^2k+2(P^Po^Q) = —Ep +PQ+q (P p+po+q Pp +po+q )(E q — E po+q ) (G9) 
for any fc. Resumming and using the identities (|199|) yields: 

P + (p,Po, 9) = — ^? + = : ]_ - P ^ P ^ ~ l (G10) 

p q lrPoq{P + Po + q) E p+po+q h+ 



P(p,pa,q) 



1 P~ P + — P~ P+ 1 R + 

L p+po+q p+q po+q q , L ^q+p 



hpoq(p + Po + q) E p+po+q E+ p q E+ 

' Epo+q - (Gil) 



Po(p + Po + q) E p+P0+q 
A simpler expression holds at coinciding points: 



/ dL dx o e- pox °- qL P OL (x o ,t\x o ,0) = = (G12) 

Jo Jo hq(po + q)E po+q E+ 

A similar calculation in the case of absorbing boundaries (AA case) gives: 

Pip, Po, q) = 7 iPq++ " Pq - P ^f^l (G13) 
lrPoq(P + Po + q)E^ +po+q E q 

as well as the semi-infinite limit L = oo with an absorbing boundary at x = 

(P + - P + P~) 
Pooo( P ,Po) = - % * +P °J o) (G14) 
lrPo{p + Po)E p+po 

At coinciding points in semi-infinite system this becomes simply 
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+oo p+ 

dx e- poXo P(x ,t\xo,0) 



In the RA case the result is 

p(p, Po , q ) = ( E t+p) p j + y q p p+ P «+ q p v+p- p q + P ^ (G15) 

p qE q lrPoq(p + q + Po)E p+po+q E- 

The AR case being obtained by the global exchange of + and — as well dbS x — y L — X 
and x — > L - x (i.e p — > -p, p — > ~Po and g — > g + p + p )- 
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